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Abstract: Alday, Gaiotto, and Tachikawa conjectured relations between certain 4d 
N = 2 supersymmetric field theories and 2d Liouville conformal field theory. We study 
generalizations of these relations to 4d theories with surface operators. For one type of 
surface operators the corresponding 2d theory is the WZW model, and for another type - 
the Liouville theory with insertions of extra degenerate fields. We show that these two 4d 
theories with surface operators exhibit an IR duality, which reflects the known relation (the 
so-called separation of variables) between the conformal blocks of the WZW model and the 
Liouville theory. Furthermore, we trace this IR duality to a brane creation construction 
relating systems of M5 and M2 branes in M-theory. Finally, we show that this duality may 
be expressed as an explicit relation between the generating functions for the changes of 
variables between natural sets of Darboux coordinates on the Hitchin moduli space. 
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1. Introduction 

One of the most interesting phenomena in supersymmetric gauge dynamics is the appearance 
of infrared (IR) duality: theories different in the ultraviolet (UV) regime may well flow to the 
same IR fixed point. A prominent example is the Seiberg duality in four-dimensional A/" = 1 
super-QCD [1]. Similar dualities exist in three dimensions [2,3] and in two dimensions [4]. 
Moreover, it is known that certain two-dimensional dualities naturally arise on the two- 
dimensional world-sheets of surface operators in four-dimensional N = 2 gauge theories [5,6]. 
In the present paper, we propose a new IR duality between 4d A/" = 2 supersymmetric theories 
with two types of surface operators that we call “codimension-2” and “codimension-4” for 
reasons that will become clear momentarily. 

In general, in four dimensional gauge theory (with any amount of supersymmetry) we 
have two ways of constructing non-local operators supported on a surface D C M4 [7]: 

• 2d-4d system: One can couple 4d gauge theory on M 4 to an auxiliary 2d theory on D 
in such a way that the gauge group G of the 4d theory is a subgroup of the global flavor 
symmetry of the 2d theory. In particular, the auxiliary 2d theory must have global 
symmetry G. 

• singularity: One replaces the four-dimensional space-time M 4 with the complement 
M 4 \ D so that gauge fields (and, possibly, other fields) have a prescribed singular 
behavior along D. Thus, instead of introducing new degrees of freedom, one modifies 
the existing degrees of freedom. 

Note that both of these methods may also be used to construct other non-local operators, 
such as line operators (for example, Wilson operators and ’t Hooft operators, respectively). 
In the case of surface operators, the first of these two methods can be further subdivided into 
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linear and non-linear sigma-model descriptions of 2d degrees of freedom on D. However, this 
distinction will not be important in this paper. 

What will be important to us, however, is that sometimes these two constructions may 
lead to the same result. This happens when integrating out 2d degrees of freedom in the 2d-4d 
coupled system leaves behind a delta-function singularity, supported on D (for the 4d fields). 
In particular, this is what one finds in the case of AA = 4 super-Yang-Mills theory. Thus, one 
obtains an equivalence of the theories with two types of surface operators, which may also be 
derived using brane constructions and T-dualities. Something similar may happen in certain 
gauge theories with less supersymmetry, e.g. free field theories, but in this paper focus on IR 
equivalence (or IR duality) of 4d AA = 2 theories with the two types of surface operators. 

Surface operators in 4d AA = 2 theories were first considered in [8] and later incorporated 
in the framework of the Alday-Gaiotto-Tachikawa (AGT) correspondence in [9,10] relating 
a certain class of 4d AA = 2 gauge theories (often called “class 5”) and 2d conformal field 
theories on a Riemann surface Cg^n of genus g with n punctures [11]. According to these 
works, there is a relation between the instanton partition functions in the 4d theories in the 
presence of the two types of surface operators and conformal blocks in the WZW model for 
SL 2 and the Liouville theory with extra degenerate fields, respectively. We note that for the 
surface operators of the first type this relation was originally proposed by Braverman [12] and 
further analyzed in [10,13-15]. 

Within this framework, the IR duality between the 4d theories with two types of surface 
operators is neatly expressed by an integral transform between the chiral partition functions 
of the WZW model and the Liouville theory: 

Z^'^{x^z) = J du K{x,u) Z^{u, z), (1-1) 

This relation, which is of interest in 2d GFT, was established by Feigin, Frenkel, and Stoy- 
anovsky in 1995 as a generalization of the Sklyanin separation of variables for the Gaudin 
model [16] (which corresponds to the limit of the infinite central charge), see [17,18]. Hence 
we call this relation separation of variables. In this paper we present it in a more explicit 
form (see [19] for another presentation). 

One of our goals is thus to show that the relation (1.1) captures the IR duality of 4d 
N = 2 gauge theories of class S with surface operators. Thus, our work provides a physical 
interpretation - and perhaps a natural home - for the separation of variables (1.1) in 4d 
gauge theory, as well as the corresponding 6d (0, 2) theory on the fivebrane world-volume in 
M-theory. 

Let’s talk about the latter in more detail. In the context of the AGT correspondence 
and, more broadly, in 4d A/" = 2 theories constructed from M-theory fivebranes wrapped on 
Riemann surfaces [20-23] the two types of surface operators in 4d field theories described 
above are usually represented by different types of branes / supersymmetric defects in the 
6d (0, 2) theory on the fivebrane world-volume. Godimension-4 defects that correspond to 
the membrane boundaries naturally lead to the surface operators described as 2d-4d coupled 


- 3 - 


systems. Codimension-2 defects, on the other hand, may be thought of as intersections with 
another group of hvebranes and therefore they are usually characterized by a singularity for 
the gauge fields at L) of a specific type (described in Appendix A). 

Thus, altogether one has at least three different perspectives on the surface operators in 
4d theories corresponding to the codimension-2 and codimension-4 defects in 6d theory (this 
is the reason why we will often refer to them as codimension-2 and codimension-4 surface 
operators). Namely, the 2d CFT perspective, the 4d gauge theory perspective, and the 
6d fivebrane / M-theory perspective. Moreover, the 4d gauge theory perspective is further 
subdivided into UV and IR regimes. A simple way to keep track of these perspectives is to 
think of a sequence of RG flows, 

M-theory / 6d 4d gauge theory UV 4d gauge theory IR (1-2) 

where arrows correspond to integrating out more and more degrees of freedom. This rela¬ 
tion between different theories is somewhat analogous to a more familiar relation between a 
2d gauged linear sigma-model, the corresponding non-linear sigma-model, and the Landau- 
Ginzburg theory that describes the IR physics of the latter. 

It is natural to ask whether one can see any trace of our IR equivalence in the UV, 
either in 4d or 6d. We answer this question in the affirmative, by showing that the brane 
configurations in M-theory that give rise to the codimension-2 and codimension-4 surface 
operators are related by a certain non-trivial phase transition, a variant of the hrane creation 
effect of Hanany and Witten [24] (see Figure 1 in Section 4.1). We will show that certain 
quantities protected by supersymmetry remain invariant under this phase transition, thereby 
revealing the 6d / M-theory origin of our IR equivalence. In four dimensions, the IR duality 
manifests itself in the most direct way as a relation between instanton partition functions in 
the presence of surface operators and conformal blocks in WZW/Liouville GFTs discussed 
above. However, what we actually claim here is that the IR duality holds for the full physical 
theories (and not just for specific observables); that is to say, the 4d theories with two types 
of surface operators become equivalent in the IR. This has many useful implications (and 
applications), far beyond a mere relation between the instanton partition functions. 

In order to show that, we use the fact that the low-energy effective action in our theories 
is essentially determined by their respective effective twisted superpotentials (see Sections 
4.2 and 4.3 for more details). Hence we need to compare the twisted superpotentials arising 
in our theories, and we compute them explicitly using the corresponding 2d conformal field 
theories. The result is that the two twisted superpotentials, which we denote by >V^®(a, x, r) 
and >V^^(a,u(a,X,r),r), respectively, are related by a field redefinition^ 

W'^®(a, x,r) = VV’^^(a, u(a, X, r), r)-I->V®°^(x, u(a, X, r), r). (1.3) 

Here the variables x and u are parameters entering the UV-definitions of the two types of 
surface operators. The relation u = u{a, x, r) extremizes the superpotential on the right of 
(1.3), reflecting the fact that u becomes a dynamical field in our brane creation transition. 

'^As usual, it is convenient to think of parameters as background fields [25]. 
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Formula (1.3) has an elegant interpretation in terms of the mathematics of the Hitchin 
integrable system for the group SL 2 . Namely, we show that the two effective twisted superpo- 
tentials are the generating functions for changes of variables between natural sets of Darboux 
coordinates for the Hitchin moduli space AiuiC) of SL 2 . 

There are in fact three such sets: {x,p), the natural coordinates on M^{C) arising from 
its realization as a cotangent bundle; (o, t), the action-angle coordinates making the complete 
integrability of AfnlC*) manifest; and {u,v), the so-called “separated variables” making the 
eigenvalue equations of the quantized Hitchin systems separate. We show that the twisted 
superpotentials x, r) and u, r) are the generating functions for the changes of 

Darboux coordinates {x,p) -H- (a, t) and {u,v) •(-)• (a, t), respectively. The generating function 
of the remaining change {x,p) -f-)- (u, v) is the function u, r) appearing on the RHS of 

the relation (1.3) - it is the generating function for the separation of variables in the Hitchin 
integrable system. 

{x,p) coordinates 



{u, v) coordinates 

Thus, the IR duality between the 4d gauge theories with the two types of surface operators 
that we study in this paper becomes directly reflected in the separation of variables of the 
Hitchin integrable system. 

To derive the relation (1.3), we first express the twisted superpotentials W^®(a,a;,T) 
and yV’'^®(a, u, r) as the subleading terms in the expansion of the logarithms of the instanton 
partition functions in the limit of vanishing Omega-deformation [26]. Assuming that the 
instanton partition function in our 4d theories are equal to the chiral partition functions in 
the WZW model and the Liouville theory, respectively [10,12-15], we express the subleading 
terms of the instanton partition functions as the subleading terms of the chiral partition 
functions in the corresponding 2d CFTs. What remains to be done then is to find a relation 
between the subleading terms of these two chiral partition functions (one from the WZW 
model and one from the Liouville theory with extra degenerate fields). 

This is now a problem in 2d CFT, which is in fact a non-trivial mathematical problem 
that is interesting on its own right. In this paper, by refining earlier observations from [27], we 
compute explicitly the subleading terms of the chiral partition functions in the WZW model 
and the Liouville theory (with extra degenerate fields) and identify them as the generating 
functions for the changes of Darboux coordinates mentioned above. In this way we obtain 
the desired relation (1.3). 
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The details of these computations are given in the Appendices, which contain a number 
of previously unpublished results that could be of independent interest. In performing these 
computations, we addressed various points in the mathematics of the WZW model and its 
relation to the Hitchin integrable system that, as far as we know, have not been discussed in 
the literature before (for example, questions concerning chiral partition functions on Riemann 
surfaces of higher genus). In particular, our results make precise the sense in which Liouville 
theory and the WZW model both appear as the result of natural quantizations of the Hitchin 
integrable systems using two different sets of Darboux coordinates, as was previously argued 
in [27]. 

Once we identify the subleading terms of the chiral partition functions of the two 2d 
CFTs with the generating functions, we obtain the relation (1.3). Alternatively, this relation 
also appears in the inhnite central charge limit from the separation of variables relation (1.1) 
between conformal blocks in the WZW and Liouville CFTs. Therefore, the relation (1.1) may 
be viewed as a relation between the instanton partition functions in the 4d theories with two 
types of surface operators in non-trivial Omega-background. This suggests that these two 
4d theories remain IR equivalent even after we turn on the Omega-deformation. However, in 
non-zero Omega-background this relation is rather non-trivial, as it involves not just a change 
of variables, but also an integral transform. This relation deserves further study, as does the 
question of generalizing our results from the group SL 2 to groups of higher rank. 


The paper is organized as follows. In Section 2 we review class S supersymmetric gauge 
theories, AGT correspondence, surface operators, and the Hitchin system. In Section 3 we 
discuss the 4d theories with the surface operators obtained from codimension-2 defects in 6d, 
the brane construction, conformal blocks in the corresponding CFT (WZW model), and the 
relation to the Hitchin system. In Section 4 we consider the 4d theories with the surface 
operators obtained from codimension-2 defects in 6d and the corresponding CFT (Liouville 
theory with degenerate fields). We also discuss general properties of the 4d theories in the 
IR regime and the corresponding twisted superpotentials. Anticipating the IR duality that 
we establish in this paper, we start with the brane system introduced in Section 3 (the one 
giving rise to the codimension-2 defects) and deform it in such a way that the end result is 
a collection of codimension-4 defects. This allows us to demonstrate that the two types of 
defects preserve the same subalgebra of the supersymmetry algebra and to set the stage for 
the IR duality. In the second half of Section 4, we bring together the results of the previous 
sections to demonstrate the IR duality of two 4d gauge theories with surface operators and 
the separation of variables in conformal field theory and Hitchin system. 

The necessary mathematical results on surface operators, on chiral partition functions in 
the WZW model and the Liouville theory, and on the separation of variables are presented 
in the Appendices. There one can also find detailed computations of the chiral partition 
functions of the WZW model and the Liouville theory and their classical limits (some of 
which have not appeared in the literature before, as far as we know). 
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2. Preliminaries 

In this section we review some background and introduce the notation that will be used in our 
paper. Toward this end, we will recall the notion of class S supersymmetric gauge theories 
and review very briefly how the Seiberg-Witten theory of this class is related to the Hitchin 
system. 

2.1 Theories of class S and AGT correspondence 

A lot of progress has been made in the last few years in the study of AA = 2 supersymmetric 
field theories in four dimensions. Highlights include exact results on the expectation values 
of observables like supersymmetric Wilson and’t Hooft loop operators on the four-sphere 5^, 
see [28,29] for reviews, and [30] for a general overview containing further references. 

A rich class of field theories with N = 2 supersymmetry, often denoted as class 5, can 
be obtained by twisted compactification of the six-dimensional (2, 0) theory with Lie algebra 
g [23]. Glass S theories of type g = Ai have Lagrangian descriptions specified by a pair of 
pants decompositions of C, which is dehned by cutting C along a system C = {71,..., 7^} of 
simple closed curves on C [22]. In order to distinguish pants decompositions that differ by 
Dehn twists, we will also introduce a trivalent graph T inside C such that each pair of pants 
contains exactly one vertex of T, and each edge e of T goes through exactly one cutting curve 
7e G C. The pair a = (C,r) will be called a refined pants decomposition. 

Then, to a Riemann surface C of genus g and n punctures one may associate [22, 23] 
a four-dimensional gauge theory Qc with N = 2 supersymmetry, gauge group (SU(2))^, 
h := 2>g — 2> + n and flavor symmetry (SU(2))”. The theories in this class are UV-finite, 
and therefore they are characterized by a collection of gauge coupling constants gi,... ,gh. 
To the /c-th boundary there corresponds a flavor group SU{2)k with mass parameter M^. 
The hypermultiplet masses are linear combinations of the parameters nik, k = l,...,n as 
explained in more detail in [ 11 , 22 ]. 

The correspondence between the data associated to the surface C and the gauge theory 
Qc is then summarized in the table below. 

We place this in the context of M-theory, following the standard conventions of brane 
constructions [20]. Namely, we choose and as local coordinates on the Riemann surface 
C and parametrize the four-dimensional space-time M 4 by (x^, x^, x^, x^). This choice of local 
coordinates can be conveniently summarized by the diagram: 



Brane 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Mb 

X 

X 

X 

X 



X 




X 


where each “x” represents a space-time dimensions spanned by the five-brane world-volume. 

Alday, Gaiotto, and Tachikawa (AGT) observed that the partition functions of Ai theories 
on a four-sphere can be expressed in terms of Liouville correlation functions. 


Riemann surface C 

Gauge theory Qc 

Cut system C -|- trivalent 

Lagrangian description with 

graph T on C, fj = (C , T) 

action functional 

cutting curve 7 e 

vector multiplet (pe, ■ ■ ■) 

n boundaries 

n hypermultiplets 

Gluing parameters qe = 

e = 1,..., h, h := 35 ^ — 3 -|- n 

UV-couplings T = {ti, ... ,Th), 

4:Tri 9e 

Te = ^ 

9 i 27r 

Change of pants decomposition 

various dualities 


2.2 Seiberg-Witten theory 

The low-energy effective actions of class S theories are determined as follows. Given a 
quadratic differential t on C one defines the Seiberg-Witten curve Ssw in T*C as follows: 

Ssw = { {u, v) G T*C ; + t{u) = 0 } . ( 2 . 1 ) 

The curve Ssw is a two-sheeted covering of C with genus 4g' — 3 -|- n. One may embed the 
Jacobian of C into the Jacobian of Ssw by pulling back the holomorphic differentials on C un¬ 
der the projection Sgw —> C. Let = LJi(Ssw,^)/J^i(G, Z), and let us introduce a 

canonical basis IB for ^)) represented by a collection of curves (ai, ..., ah; q;°, ..., a^) 

with intersection index o af = 5 ^ 1 , a^oai = 0, a^ o af = 0. The corresponding periods of 
the canonical differential on u = v{u)du are defined as 

ak = f V, = f V. (2.2) 

Using the Riemann bilinear relations, it can be shown that there exists a function T'(a), 
a = (oi,..., ah) such that = da,.F{a). The function F{a) is the prepotential determining 
the low-energy effective action associated to B. 

Different canonical bases B for LJ((Ssw)^) are related by S'p(2h, Z)-transformations de¬ 
scribing electric-magnetic dualities in the low-energy physics. It will be useful to note that for 
given data a specifying UV-actions there exists a preferred class of bases B^ for iL((Ssw,Z) 
which are such that the curves project to the curves 7 eGC,e = l,...,/i defining the pants 
decomposition C, respectively. 



2.3 Relation to the Hitchin system 

The Seiberg-Witten analysis of the theories Qc has a well-known relation to the mathematics 
of the Hitchin system [31,32] that we will recall next. 

The phase space of the Hitchin system for G = SL{2) is the moduli space of 

pairs where T is a holomorphic rank 2 vector bundle with fixed determinant, and 

If £ H^{C, End(T) ( 8 > Kq) is called the Higgs field. The complete integrability of the Hitchin 
system is demonstrated using the so-called Hitchin map. Given a pair {£,ip), we define the 
spectral curve S as 

S = { (u, u) G T*C ; 2^2 = ti{ip^{u)) } . (2.3) 

To each pair (T, ip) one associates a line bundle L on S, the bundle of eigenlines of f for a 
given eigenvalue v. Conversely, given a pair (S,L), where S C T*C is a double cover of C, 
and L a holomorphic line bundle on S, one can recover (£, ip) via 

{S,ip) := (7r*(L), 7r*(u)) , (2.4) 

where tt is the covering map S —)■ C, and tt* is the direct image. 

The spectral curves S can be identified with the curves Ssw determining the low-energy 
physics of the theories Gc on However, in order to give physical meaning to the full 
Hitchin system one needs to consider an extended set-up. One possibility is to introduce 
surface operators. 

2.4 Two types of surface operators 

When the 6 d fivebrane world-volume is of the form M 4 x C, where (7 is a Riemann surface, 
there are two natural ways to construct half-BPS surface operators in the four-dimensional 
space-time M 4 where the N = 2 theory Qc lives. First, one can consider codimension-2 
defects supported on D x C, where D C M4 is a two-dimensional surface (= support of a 
surface operator). Another, seemingly different way, is to start with codimension-4 defects 
supported on D x {p}, where p £ C is a point on the Riemann surface. 

In the case of genus-1 Riemann surface C = T‘^, both types of half-BPS surface operators 
that we study in this paper were originally constructed using branes in [7,33]. In these 
papers it was argued that the two types of operators are equivalent, at least for certain 
“supersymmetric questions”. Here we will show that for more general Riemann surfaces 
C the two surface operators, based on codimension-4 and codimension -2 defects, may be 
different in the UV but become essentially the same in the IR regime. They correspond 
to two different ways to describe the same physical object. Mathematically, this duality of 
descriptions corresponds to the possibility of choosing different coordinates on the Hitchin 
moduli space, which will be introduced shortly. At first, the equivalence of the two types 
of surface operators may seem rather surprising since it is not even clear from the outset 
that they preserve the same subalgebra of the supersymmetry algebra. Moreover, the moduli 
spaces parametrizing these surface operators appear to be different. 
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Indeed, one of these moduli spaces parametrizes collections of n codimension-4 defects 
supported at D x {pi} C M 4 x C, and therefore it is 

Sym’^(C') := /Sn (2.5) 

(Here we consider only the “intrinsic” parameters of the surface operator, and not the position 
of D C. M4, which is assumed to be fixed.) On the other hand, a surface operator constructed 
from a codimension -2 defect clearly does not depend on these parameters, since it wraps on 
all of C. Instead, a codimension -2 surface operator carries a global symmetry G — which 
plays an important role e.g. in describing charged matter — and, as a result, its moduli space 
is the moduli of G-bundles on C, 

BunG(C) (2.6) 

Therefore, it appears that in order to relate the two constructions of surface operators, one 
must have a map between (2.5) and (2.6): 

BunG(C') ^ Sym"(C) ^^. 7 ) 

X ^ u 

where n = {g — l)dimG = dimBunG'(C'). 

It turns out that even though such a map does not exist, for G = SL{2) there is a map 
of the corresponding cotangent bundles, which is sufficient for our purposes. This is the 
celebrated classical separation of variables. Moreover, it has a quantum version, described in 
Section 4.7. The separation of variables allows us to identify the 4d theories with two types 
of surface operators in the IR. 

The unbroken SUSY makes it possible to turn on an Omega-deformation, allowing us to 
define generalizations of the instanton partition functions. In the case of codimension-2 surface 
operators it turned out that the generalized instanton partition functions are calculable by 
the localization method, and in a few simple cases it was observed that the results are related 
to the conformal blocks in the SL{2)-WZW model. For codimension-4 surface operators 
one expects to find a similar relation to Liouville conformal blocks with a certain number of 
degenerate fields inserted. 

3. Surface operators corresponding to the codimension-2 defects 

Our goal in this paper is to establish a relation between the surface operators constructed 
from codimension-2 and codimension-4 defects.^ In order to do that, we must show that they 
preserve the same subalgebra of the supersymmetry algebra. This will be achieved by realizing 
these defects using branes in M-theory (as we already mentioned earlier). This realization 
will enable us to link the two types of defects, and it will also illuminate their features. 

^Even though our main examples will be theories of class S, we expected our results — in particular, the 
IR duality — to hold more generally. 
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In this section we present an M-theory brane construction of the codimension-2 defects 
and then discuss them from the point of view of the 4d and 2d theories. Then, in Section 4, 
we will deform — in a way that manifestly preserves supersymmetry — a brane system that 
gives rise to the codimension-2 defects into a brane system that gives rise to codimension-4 
defects. Using this deformation, we will show that the two types of defects indeed preserve 
the same supersymmetry algebra, and furthermore, we will connect the two types of defects, 
and the corresponding 4d surface operators, to each other. 

3.1 Brane construction 

Following [7], we denote the support (resp. the fiber of the normal bundle) of the surface 
operator inside M 4 by D (resp. D'). In fact, for the purposes of this section, we simply take 
M 4 = D X D'. Our starting point is the following “brane construction” of 4d AA = 2 gauge 
theory with a half-BPS surface operator supported on D C {= D x D'): 

M5 : DxD' xC 

M5' : DxC X D" (3.1) 

embedded in the eleven-dimensional space-time DxD'xT*CxMx D" in a natural way. 
For simplicity, we will assume that D = D' = D” = and C is the only topologically non¬ 
trivial Riemann surface in the problem at hand. And, following the standard conventions 
of brane constructions [ 20 ], we use the following local coordinates on various factors of the 
eleven-dimensional space-time: 


D 

D' 

T*C 

M 

D" 


2 

X ^ tX ^ X ^ X 


8 9 


With these conventions, the brane configuration (3.1) may be equivalently summarized in the 
following diagram: 


Brane 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

M5 

X 

X 

X 

X 



X 




X 

M5' 

X 

X 





X 


X 

X 

X 


Note that M5'-branes wrap the same UV curve C as the M5-branes. This brane configuration 
is |-BPS, i.e. it preserves four real supercharges out of 32. Namely, the eleven-dimensional 
space-time (without any fivebranes) breaks half of supersymmetry (since T*C is a manifold 
with SU{2) holonomy), and then each set of fivebranes breaks it further by a half. 

In particular, thinking of T*C as a non-compact Calabi-Yau 2-fold makes it clear that 
certain aspects of the system (3.1), such as the subalgebra of the supersymmetry algebra 
preserved by this system, are not sensitive to the details of the support of M5 and M5' branes 
within T*C as long as both are special Lagrangian with respect to the same Kahler form uj 
and the holomorphic 2-form U. Since T*C is hyper-Kahler, it comes equipped with a sphere 
worth of complex structures, which are linear combinations of I, J, K, and the corresponding 
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Kahler forms uij, ujj, cok- Without loss of generality, we can choose u! = uj and = ujj+iuK- 
Then, the special Lagrangian condition means that both ojj and ojk vanish when restricted 
to the world-volume of M5 and M5' branes. 

3.2 Four-dimensional description 

As we explain below, surface operators originating from codimension-4 defects in 6 d (0, 2 ) 
theory naturally lead to the coupled 2d-4d system, while those originating from codimension- 
2 defects in 6 d descend to the second description of surface operators in 4d gauge theory, 
namely as singularities for the UV gauge fields (see Appendix A for more details); 

AirW- (*)’■' (3.3) 

Here, following our conventions (3.2), we use a local complex coordinate + ix^ = r 2 e*®^ on 
D' such that surface operator is located at the origin (r 2 = 0). A surface operator defined 
this way breaks half of supersymmetry and also breaks SO {A) rotation symmetry down to 
50(2) X 50(2). From the viewpoint of the 2d theory on D, the unbroken supersymmetry is 
AA= (2,2). 

The symmetries preserved by such a surface operator are exactly what one needs in order 
to put the 4d gauge theory in a non-trivial Omega-background. Mathematically, this leads 
to an 50(2) x 50(2) equivariant counting of instantons with a ramification along D. The 
resulting instanton partition function 

Z’^®(a,x,r;ei,e2), (3.4) 

depends on variables x = (xi,... ,Xh) related to the parameters in (3.3) via the expo¬ 
nentiation map 

Xr = . (3.5) 

The relation between the parameters and the counting parameters Xr appearing in the 
instanton partition functions was found in [ 10 ]. 

3.3 Relation to conformal field theory 

Starting from the groundbreaking work of A. Braverman [12], a number of recent studies have 
produced evidence of relations between instanton partition functions in the presence of surface 
operators Z^^{a, x, r; ei, € 2 ) and conformal blocks of affine Kac-Moody algebras Qk [10,13-15]. 
Such relations can be viewed as natural generalizations of the AGT correspondence. In the 
case of class 5-theories of type Ai one needs to choose g = 5 I 2 and k = —2 — as will be 
assumed in what follows. 

The Lie algebra Qk has generators J“, a = 0, , n G Z. A large class of representation 

of Qk is defined by starting from a representation Rj of the zero mode subalgebra generated 
from Jq, which has Casimir eigenvalue parametrized as j{j -|- 1). One may then construct a 
representation TZj of g^ as the representation induced from TZj extended to the Lie subalgebra 
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generated by J“,n > 0, such that all vectors v € Rj C TZj satisfy J^v = 0 for n > 0. To be 
specific, we shall mostly discuss in the following the case that the representations TZj have a 
lowest weight vector Cj, but more general representations may also be considered, and may 
be of interest in this context [34], 

In order to define the space of conformal blocks, let (7 be a compact Riemann surface and 
zi,..., an n-tuple of points of C with local coordinates ti,..., We attach representations 

TZr = TZj^ of the affine Kac-Moody algebra gfc of level k to the points Zr-, r = ... ^n. The 

diagonal central extension of the direct sum 0"=^ 0 ® C((tr)) acts on the tensor product 
Consider the Lie algebra 

0 out = 0 (8) C[C'\{zi, . . . , Zn}] 

of 0-valued meromorphic functions on C with poles allowed only at the points zi,..., We 
have an embedding 

n 

0 out ^ ( 3 - 6 ) 

r=l 

It follows from the commutation relations in 0 and the residue theorem that this embed¬ 
ding lifts to the diagonal central extension of 0”=;^ 0 0 C{{tr))- Hence the Lie algebra 0out 
acts on (S)r=i'^r- By definition, the corresponding space of conformal blocks is the space 
CBg(7?.i,... ,TZn) of linear functionals 

n 

(f : TZ^^n] ■= TZr ^ C 

r=\ 

invariant under 0 out, i-e., such that 

n 

if{ri-v) = 0, yvel^TZr, ri e Q®C[C\{zi,...,Zn}]. (3.7) 

r=l 

The conditions (3.7) represent a reformulation of current algebra Ward identities well-known 
in the physics literature. The space CBg(77i,... ,77„) is infinite-dimensional in general. 

To each ip S CBg(77i,..., TZn) we may associate a chiral partition function Z{ip, C) by 
evaluating ip on the product of the lowest weight vectors, 

Z^'^{p,C]k) \= Cn). (3.8) 

In the physics literature one usually identifies the chiral partition functions with expectation 
values of chiral primary fields ^r{zr), inserted at the points 

C-k) = { ^n{Zn) • • • 4>l(^l) . (3.9) 

Considering families of Riemann surfaces 0 parametrized by local coordinates r for the 
Teichmliller space Tg^n one may regard the chiral partition functions as functions of r, 

Z'^\p,Cpk) = Z^\p,T-k). 
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Large families of conformal blocks and the corresponding chiral partition functions can 
be constructed by the gluing construction. Given a (possibly disconnected) Riemann surface 
C with two marked points Pq , i = 1,2 surrounded by parametrized discs Dj one can construct 
a new Riemann surface by pairwise identifying the points in annuli Aj C Oj around the two 
marked points, respectively. Assume we are given conformal blocks associated to two 
surfaces C* with n* + 1 punctures Pq, ..., P^. with the same representation TZq associated 
to Pq for i = 1,2. Using this input one may construct a conformal block <^<7^2 associated to 
the surface C 12 obtained by gluing the annular neighborhoods Aj of Pg, i = 1, 2 as follows: 

Aci2 (^^1 ® • • • 'S’Vni (8) rci <8) • • • <8) Wn2 ) = 

= ^ 'Si v,y) iPq^{)C{t,x)v'^'Si Wi'Si ■ ■ ■'Si Wn2) ■ 

The vectors and are elements of bases for the representation TZq which are dual w.r.t. to 
the invariant bilinear form on TZq. A standard choice for the twist element /C(t, x) G End(Po) 
appearing in this construction is /C(t, x) = , where the operator Lq represents the 

zero mode of the energy-momentum tensor constructed from the generators J“ using the 
Sugawara construction. The parameter q = in (3.10) can be identified with the modulus 
of the annular regions used in the gluing construction of Ci 2 - However, it is possible to 
consider twist elements IC{t,x) constructed out a larger subset of the generators of gfc. The 
rest of the notation in (3.10) is self-explanatory. The case that Pg, i = 1,2 are on a connected 
surface can be treated in a similar way. 

A general Riemann surface Cg^n can be obtained by gluing 2g — 2+ n pairs of pants Cq 3, 
V = 1,... ,2g — 2 + n. It is possible to construct conformal blocks for the resulting Riemann 
surface from the conformal blocks associated to the pairs of pants Cg 3 by recursive use of the 
gluing construction outlined above. This yields families of conformal blocks parametrized 
by 

• the choice of a refined pants decomposition a = (C,r), 

• the choice of representation TZj^ for each of the cutting curves 7e defined by the pants 
decomposition, and 

• the collection of the parameters Xg introduced via (3.10) for each curve 7e G C. 

The corresponding chiral partition functions are therefore functions 

X, r; k) = r; k). 

The variables x = (xi,... ,X 3g_3+„) have a geometric interpretation as parameters for fam¬ 
ilies of holomorphic G = S'P(2)-bundles B. Indeed, in Appendix B it is explained how the 
definition of the conformal blocks can be modified in a way that depends on the choice of 
a holomorphic bundle B, and why the effect of this modihcation can be described using the 
twist elements /C(r, x) appearing in the gluing construction. It follows from the discussion in 
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Appendix B that changing the twist elements /C(t, x) amounts to a change of local coordinates 
(r, x) for the fibration of Bunc over Tg^n (the moduli space of pairs: a Riemann surface and 
a G-bundle on it). 

The chiral partition functions satisfy the Knizhnik-Zamolodchikov-Bernard (KZB) equa¬ 
tions. This is a system of partial differential equations of the form 


€2 d 
ei dqe 


zrij, 


X, r; k) 


^e2^^U,x,T-,k), 


(3.11) 


where Hg is a second order differential operator containing only derivatives with respect to 
the variables Xg- These equations can be used to generate the expansion of x; r; k) in 

powers of qg and Xg, 


Za^{j,x,T]k) ~ 


E 

nezi 


27^U,m,n-,k)Y{q 




(3.12) 


mez5 


e=l 


The notation ~ used in (3.12) indicates equality up to a factor which is j-independent. Such 
factors will be not be of interest for us. The equations (3.11) determine Z^^(j,in,n;k) 
uniquely in terms of Z^^{j) = Z'^'^{j, 0, 0; k). It is natural to assume that the normalization 
factor Z^^{j) can be represented as product over factors depending on the choices of repre¬ 
sentations associated to the three-holed spheres Cq 3 appearing in the pants decomposition. 

We are now going to propose the following conjecture: There exists a choice of twist 
elements JCg{Tg,Xg) such that we have 


Zf\a,x,T-ei,€ 2 ) ^ Zj^{j,x,q-,k), 


(3.13) 


assuming that 

je = -^+i—, k + 2 = -—. (3.14) 

2 Cl €1 

Evidence for this conjecture is provided by the computations performed in [10,13-15] in the 
cases C = Ciq and C = The relevant twist elements JC{t,x) were determined explicitly 
in these references. As indicated by the notation we expect (3.13) to hold only up to 
j-independent multiplicative factors. A change of the renormalization scheme used to define 
the gauge theory under consideration may modify Z^^ by factors that do not depend on j. 
Such factors are physically irrelevant, see e.g. [35] for a discussion. 


3.4 Relation to the Hitchin system 

On physical grounds we expect that the instanton partition functions Z^^{a,x,T]ei,€ 2 ) be¬ 
have in the limit ei —)> 0, €2 —> 0 as 

log2:^®(a,x,T; € 1 , 62 ) ~- —Ta{a,T)-—W^^{a,x,T). (3.15) 

£162 £1 

The first term is the bulk free energy, proportional to the prepotential defined previ¬ 

ously. The second term is a contribution diverging with the area of the plane on which the 
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surface operator is localized. It can be identified as the effective twisted superpotential of the 
degrees of freedom localized on the surface X2 = X3 = 0. 

The expression of the instanton partition function as a to conformal field theory (3.13) 
allows us to demonstrate that we indeed have an asymptotic behavior of the form (3.15). The 
derivation of (3.15) described in Appendix D leads to a precise mathematical description of 
the functions W^®(o, x, r) appearing in (3.15) in terms the Hitchin integrable system that we 
will describe in the rest of this subsection. It turns ont that W’^®(a, x, r) can be characterized 
as the generating fnnction for the change of variables between two sets of Darboux coordinates 
for AIh(C') naturally adapted to the description in terms of Higgs pairs (T, (p) and pairs (S, L), 
respectively. 

Let ns pick coordinates x = (xi,..., Xh) for Bunc. Possible ways of doing this are briefly 
described in Appendix C.2. One can always hnd coordinates p on AIh(C') which supplement 
the coordinates x to a system of Darboux coordinates (x,p) for 

There exists other natural systems (a, t) of coordinates for Ai^{C) called action-angle 
coordinates making the complete integrability of AdH(C') manifest. The coordinates a = 
(ai,..., o/i) are defined as periods of the Seiberg-Witten differential, as described previously. 
The coordinates t = (ti,..., th) are complex coordinates for the Jacobian of S parametrizing 
the choices of line bundles L on S. The coordinates t may be chosen such that (a, t) furnishes 
a system of Darboux coordinates for AJh(C')- 

As the coordinates (a, t) are natnrally associated to the description in terms of pairs 
(S, L), one may construct the change of coordinates between the sets of Darboux coordinates 
{x,p) and (a, t) nsing Hitchin’s map introduced in Section 2.3. The function >V^®(a,x,T) 
in (3.15) can then be characterized as the generating function for the change of coordinates 
{x,p) ^ {a,t), 






27r dar 




(3.16) 


with periods a defined using a basis corresponding to the pants decomposition a used 
to define .Z^®(a, x, r; ei, 62). Having defined (x,p) and (a, t), the equations (3.16) define 
W^®(a,x,r) up to an (inessential) additive constant. 


3.5 Physical interpretation 

All of the integrable system gadgets introduced above seem to find natural homes in field 
theory and string theory. In particular, N five-branes on C describe a theory that in the IR 
corresponds to an M5-brane wrapped N times on C or, equivalently, wrapped on a A^-fold 
cover S —)• C. 

Though in this paper we mostly consider the case N = 2 (hence a double cover S —)• C), 
certain aspects have straightforward generalization to higher ranks. It is also worth noting 
that we treat both SL{N) and GL{N) cases in parallel; the difference between the two is 
accounted for by the “center-of-mass” tensor multiplet in 6d (0, 2) theory on the five-brane 
world-volume. 
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Besides the “brane constructions” used in most of this paper, the physics of 4d AA = 2 
theories can be also described by compactification of type IIA or type IIB string theory on a 
local Calabi-Yau 3-fold geometry. This approach, known as “geometric engineering” [36,37], 
can be especially useful for understanding certain aspects of surface operators and is related to 
the brane construction by a sequence of various dualities. Thus, a single five-brane wrapped 
on E C that describes the IR physics of 4d A/" = 2 theory is dual to type IIB string 
theory on a local CY 3-fold 

zw — P{u,v) = 0, (3-17) 

where P{u,v) is the polynomial that defines the Seiberg-Witten curve Esw 

It can be obtained from our original M5-brane on S by first reducing on one of the 
dimensions transversal to the five-brane (down to type IIA string theory with NS5-brane on 
E) and then performing T-duality along one of the dimensions transversal to the NS5-brane. 
The latter is known to turn NS5-branes to pure geometry, and supersymmetry and a few 
other considerations quickly tell us that type IIB background has to be of the form (3.17). 

Now, let us incorporate M5'-brane which in the IR version of brane configuration (3.1) 
looks like: 


M5 : D X D' X T, 

MS' : D X E X D" (3.18) 

What becomes of the MS'-brane upon duality to type IIB setup (3.17)? 

It can become any brane of type IIB string theory supported on a holomorphic subman¬ 
ifold in the local Calabi-Yau geometry (3.17). Indeed, since the chain of dualities from M- 
theory to type IIB does not touch the four dimensions parametrized by ,..., the resulting 
type IIB configuration should still describe a half-BPS surface operator in 4d Seiberg-Witten 
theory on M4. Moreover, since type IIB string theory contains half-BPS p-branes for odd 
values of p, with (p -|- l)-dimensional world-volume, MS' can become a p-brane supported on 
D X Cp-i, where Cp-i is a holomorphic submanifold in a local Calabi-Yau 3-fold (3.17). 

Depending on how one performs the reduction from M-theory to type IIA string theory 
and then T-duality to type IIB, one finds different p-brane duals of the MS'-brane. Here, we 
will be mostly interested in the case p = 3, which corresponds to the reduction and then T- 
duality along the coordinates x® and x®, cf. (3.2). Effectively, one can think of compactifying 
the M-theory setup (3.18) on D” = T^, and that gives precisely the type IIB setup (3.17) 
with extra D3-brane supported on E, i.e. aX z = w = Q m. (3.17). 

A D3-brane carries a rank-1 Chan-Paton bundle L' —)• E. Therefore, we conclude that the 
surface operators made from codimension-2 defects that are obtained from the intersections 
with MS'-branes as described above, have an equivalent description in dual type IIB string 
theory in terms of pairs (E,L'). It seems likely that the line bundle L' is closely related to 
the line bundle L appearing in the description of the Hitchin system in terms of pairs (E, L). 

Note, the degree of this line bundle, d{L'), is equal to the induced Dl-brane charge along 
the (x°,x^) directions. For completeness, we describe what it corresponds to in the dual M- 
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theory setup (3.18). The T-duality that relates type IIA and type IIB brane configurations 
maps Dl-branes supported on into D2-branes with world-volume along 

Hence, we conclude 

d{L') = M2-brane charge along X®) (3.19) 

It seems worthwhile investigate the description of surface operators in terms of type IIB brane 
configurations in more detail. 

4. Surface operators corresponding to codimension-4 defects 

As we mentioned earlier, there is another way to construct surface operators in 4d A/" = 2 
theories of class S - namely, by introducing codimension-4 defects in 6d five-brane theory 
[20-22,38]. 

In this section we present this construction. The idea is to start with the brane system 
which we used in the previous section to produce the codimension-2 defects and to deform it 
in such a way that the end result is a collection of codimension-4 defects. The advantage of 
this way of constructing them is that, as we will see below, this process does not change the 
subalgebra of the supersymmetry algebra preserved by the defects. Therefore, it follows that 
the two types of defects in fact preserve the same subalgebra. 

In the next sections we will also use this link between the codimension-4 and codimension- 
2 defects in the 6d theory in order to establish the connection between the corresponding 4d 
N = 2 theories in the IR. 

4.1 Brane construction 

The origin of codimension-4 defects in 6d theory and the resulting surface operators in 4d 
N = 2 theory are best understood via the following brane construction: 


Brane 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

M5 

X 

X 

X 

X 



X 




X 

M2 

X 

X 






X 





where in addition to N M5-branes supported on M4 x C (as in Section 3.1) we have added 
a number of M2-branes supported on H x M_|_, where M+ = {x^ > 0}. Note that each of 
these M2-branes is localized at one point of the UV curve C and therefore gives rise to a 
codimension-4 defect in the 6d theory. 

One of the main goals of this paper is to show that the surface operators in 4d A/" = 2 the¬ 
ory corresponding to these codimension-4 defects describe in the IR the same physical object 
as (3.1), up to a field transformation (which is related to a change of Darboux-coordinates in 
the associated integrable system). For such an equivalence to make sense, it is necessary that 
the two types of defects preserve the same supersymmetry subalgebra. This is a non-trivial 
statement that we explain presently. 

A simple and elegant way to analyze supersymmetry and to gain further insight into the 
relation between the two types of surface operators is to perform a continuous deformation 
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a) b) c) 

Figure 1: An M5'-brane wrapped on the curve C can be perturbed to a curve C which meets C at 
finitely many points Ui. Then, separating the five-branes on C and C along the x'^ direction results in 
creation of M2-branes (shown in red). 

of one brane configuration into the other preserving the corresponding subalgebra of the 
supersymmetry algebra.^ Starting with our original system (3.1), we keep the world-volume 
of the M5-branes to be H x x C, but deform the support of the Mbhbranes to be H x C x D”, 
where C <Z T*C is a deformation of the zero section C C T*C, which is special Lagrangian 
with respect to w = w/ and Vl = ojj + iojx '■ 

M5 : Dx D' xC 

M5' : D X C X D" (4.1) 

According to the discussion in Section 3.1, this deformation does not affect the amount of 
unbroken supersymmetry, and so (4.1) preserves the same part of the supersymmetry algebra 
as the original system (3.1). Note that deformations of special Lagrangian submanifolds are 
infinitesimally parametrized by H^{C) and, in most cases of interest, this is a fairly large 
space. However, what’s even more important is that, after the deformation, C meets the 
original curve C only at finitely many points Ui, as illustrated on Figure lb. The number of 
such intersection points is determined by the Euler characteristic (or genus) of the curve C 

C-C = 2g{C)-2. (4.2) 

At low energies one may effectively represent the stack of M5-branes in terms of a smooth 
curve Tj <Z T*C [20]. The MS'-branes will be represented by a curve S' related to S by 
holomorphic deformation. Using the same arguments as above one may show, first of all, 
that two types of IR surface operators preserve the same SUSY and, furthermore, determines 
the number of intersection points on S to be 


S • S = 2fe - 2 , (4.3) 

^The argument presented below applies equally well to a system where the UV curve C is replaced by the 
IR curve E. In fact, the latter version, which similarly explains that IR surface operators preserve the same 
SUSY is also responsible for the IR duality that underlies the separation of variables map. 
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where gj, = 4g' — 3 if C has no punctures [31], as will be assumed in this section for simplicity. 

After the deformation, every intersection of M5 and M5' locally looks like a product 
of with a submanifold in M®, which is a union of two perpendicular 4-spaces U M^, 
intersecting at one point, times the real line M parametrized by the coordinate . Indeed, 
M5 and M5' overlap along a 2-dimensional part of their world-volume, D, and the remaining 
4-dimensional parts of their world-volume span = {x^ = 0}- If we separate these five- 
branes in the direction, they become linked in the 9-dimensional space which is the part 
of the space-time orthogonal to D. Then, if we make one of the five-branes pass through the 
other by changing the value of its position in the direction, an M2-brane is created, as 
shown on Figure Ic. The support of the M2-brane is D x I, where / is the interval along x^ 
connecting the deformations of the 4-spaces, which we denote by and (where a and b 
are the values of the coordinate x^ corresponding to these two subspaces): 

M5 : DxRI 

M5' : Dx 

M2 : Dx{a<x'^ <b} (4.4) 

This creation of the M2-brane between two linked M5-branes is a variant of the so-called 
Hanany-Witten effect [24]. What this means for us is that a surface operator represented 
by a codimension-2 defect wrapped on H x S in the fivebrane theory can be equivalently 
represented by a collection of codimension-4 defects supported at various points ttj G S. 

Indeed, globally, after separating M5 and M5^ in the x^ direction, the brane configuration 
(4.1) looks like this: 

M5 : D X D' xT, 

MS' : Dx D" X S 

M2 : Dxl (4.5) 

Here, adding M2-branes does not break supersymmetry any further, so that (4.5) is a |-BPS 
configuration for arbitrary special Lagrangian submanifolds S and S C r*S. Of course, the 
special case S = S takes us back to the original configuration (3.1), schematically shown in 
Figure la. On the other hand, separating M5 and MS' farther and farther apart, we basically 
end up with the standard brane configuration, shown on Figure 26, that describes half-BPS 
surface operator(s) built from codimension-4 defects, or M2-branes. In fact, even our choice 
of space-time conventions (3.2) agrees with the standard notations used in the literature, so 
that (4.5) can be viewed as M-theory lift of the following brane system in type IIA string 
theory: 


NS5 : 

012345 

D4 : 

0123 

NS5' : 

01 45 

D2 : 

01 


(4.6) 
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Conversely, reduction of (4.5) on the M-theory circle (parametrized by gives the type 
IIA system (4.6) shown on Figure 2a. 



Figure 2: The brane construction of a surface operator in pure Af = 2 super Yang-Mills theory (a) 
in type IIA string theory and (b) its M-theory lift. 


How many M2-branes are created in the configuration (4.5)? If the number of M5-branes 
is N and the number of M5'-branes is k, then each intersection point rt* G S n S contributes 
k ■ N M2-branes (due to the s-rule [24]). When we multiply this by the number of intersection 
points (4.3), we get the answer 2{g — l)kN. This number, however, counts how many M2- 
branes are created as one pulls a stack of M5'-branes through the stack of M5-branes by 
changing their x^-position from x'^ < 0 to x^ > 0 , while we are interested in a process that 
starts at x^ = 0 and then goes to either x^ < 0 or x^ > 0 . 

The initial value x^ = 0 is somewhat singular. However, as in a similar “geometric 
engineering” of 2d field theories with the same amount of supersymmetry [39], we shall assume 
that both phases x^ < 0 and x^ > 0 are symmetric and the same number of M2-branes is 
created (or destroyed) as we pass from x^ = 0 to either x^ < 0 or x^ > 0. In fact, via a chain 
of dualities [40] our “brane engineering” of the 2d theory on M2-branes can be mapped to 
the “geometric engineering” of [39], which therefore justifies applying the same arguments. 
Then, it means that the answer we are looking for is only half of 2{g — 1)A:A^, i.e. 

#(M2-branes) = {g — l)kN (4.7) 

The case considered in this paper is = A: = 2, giving a number of 4g — 4 M2-branes created. 

In the IR one may represent the M5 by a curve T, in T*C. The M5'-branes are supported 
on a holomorphic deformation of S, which may be represented by a section of a line bundle 
of the same degree as ATs, 

deg(As) = 2g^-2 = 2{Ag - 3) - 2 = 85 - 8 . (4.8) 
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It seems natural to assume that S' is symmetric under the involution exchanging the two 
sheets of S. This implies that the projection vr : S —)■ C of the intersection points defines 
Ag — A points u = (ui,... on C. Following the discussion above, one expects to find 

a collection of M2-branes created with end-points at Ur, r = 1, ... ,4g — A. 

Since a surface operator supported on C M4 breaks translation invariance in the 
transverse directions (along D'), it must necessarily break at least part of supersymmetry 
of the Ad N = 2 gauge theory on M4. In addition, our analysis above shows that both 
types of surface operators preserve the same part of supersymmetry. It is convenient to 
express the unbroken parts of 4d Lorentz symmetry and supersymmetry in 2d language. 
Indeed, the unbroken generators of the Lorentz symmetry (in and directions along 
D) conveniently combine with the unbroken supercharges and the R-symmetry generators to 
form 2 d N = (2, 2) supersymmetry algebra. 

4.2 Four-dimensional description 

We now start discussing the implications of this construction for the IR physics of 4d AA = 2 
gauge theories with surface operators. 

The Lagrangian of a 4d W = 2 gauge theory with surface operators may have additional 
terms corresponding to 2d W = (2, 2) supersymmetric theories coupled to the surface opera¬ 
tors. Recall that the Lagrangian of a theory with 2d N = (2, 2) supersymmetry is allowed to 
have a particular type of F-term called the twisted superpotential, denoted by W. From the 
point of view of a 4d theory, such a term is a two-dimensional feature, i.e. such terms would 
not be present in a 4d AA = 2 theory without surface operators, and it is partially protected 
by the supersymmetry from quantum corrections. Moreover, in the IR, the Ad N = 2 gauge 
theory with surface operators is completely determined by the prepotential J- and the twisted 
superpotential W (see e.g. [41] for a recent review). 

Recall that the low-energy effective action has a four-dimensional part and a two-dimen¬ 
sional part, 

S = j d^xd^9T+ Q j dW^W + c.c.^ , (4.9) 

where J- is the prepotential giving the low-energy effective action of the four-dimensional 
theory in the absence of a surface operator, and W is the holomorphic twisted superpotential. 
We will mostly consider A" as a function A(a,r), with a being a collection a = (ai,... ,ah) 
of coordinates for the moduli space of vacua Advac, where h is the dimension of Ad vac, and 
r being the collection of UV gauge coupling constants r = (n,..., r/i). The dependence on 
the mass parameters will not be made explicit in our notations. W = W{a,K,T) depends 
on a and r, and may furthermore depend on a collection of parameters k characterizing the 
surface operator in the UV. 

The presence of surface operators implies that the abelian gauge fields Ar, r = 1,... ,h 
appearing in the same vector-multiplet as the scalars will generically be singular at the 
support D of the surface operator. The singularity is such that the field strength associated 
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to Ar has a singularity of the form {Fr )23 = 2TrarS{x2)S{x3). The parameters Or are related 
to the twisted superpotential W by a relation of the form 


ir — T .— 


1 d 

27r dttr 


W, 


'^rs 


d d 

dar dag 


T. 


(4.10) 


The parameters r/r in (4.10) characterize the divergence of the dual gauge fields in a similar 
way. As indicated in (4.10), it is useful to combine the Gukov-Witten parameters ar and rjr 
into complex variables t = {ti,... ,th) which are functions of a, r and n. 

The argument of the previous subsection shows that the brane configuration (3.1) that 
describes codimension-2 defects can be continuously deformed without changing the unbroken 
supersymmetry to a brane configuration describing codimension-4 defects: 


M5 : Dx D' xC 

M2: DxR+ (4.11) 


This has important implications for our story. First, it means that the same type of Omega- 
background in both cases leads to the same kind of F-terms (appearing in the instanton 
partition functions) for both types of surface operators. Namely, in the language of unbroken 
2d J\f = (2,2) supersymmetry, it is the twisted superpotential W in both (3.15) and (4.18). 

Note that by itself, the existence of a continuous deformation relating surface operators 
corresponding to the codimension-2 defects to those corresponding to the codimension-4 de¬ 
fects does not necessarily imply their equivalence. Indeed, there are many physical systems 
related by a continuous deformation which describe completely different physics, e.g. gauge 
theory at different values of a coupling constant is a simple example. However, certain quan¬ 
tities may be insensitive to a change of parameter, and in fact, in the case at hand, we will 
show that the twisted superpotential W is precisely such a quantity that does not depend on 
the deformation described in the previous subsection (up to a change of variables). 

But the twisted superpotential W determines the vacuum structure and the IR physics 
of the 4d theories with surface operators. Therefore if we can show that W is independent of 
the deformation, it will follow that the corresponding 4d theories are equivalent in the IR. 

So, our plan is the following. In this subsection, we show that the twisted superpotential 
W is indeed independent of the separation of M5 and M5' in the direction, which was our 
deformation parameter in the brane configuration (4.5) that interpolates between (3.1) and 
(4.11). And then, in the next section, we will use this independence of W on the deformation 
parameter to argue that the 4d theories with the surface operators corresponding to the 
codimension-2 and codimension-4 defects describe the same physics in the IR regime (in 
other words, they are related by an IR duality). 

In order to show the x^-independence of W, we need to focus more closely on the surface 
operators produced from codimension-4 defects and explain a few facts about the brane 
systems (4.5)-(4.11) that involve M2-branes. As we already pointed out earlier, the brane 
configuration (4.5) is simply an M-theory lift of the brane system (4.6) illustrated in Figure 
2a. Usually, such M-theory lifts capture IR quantum physics of the original type IIA system, 
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cf. [20]. In the present case, the relevant theory “lives” on D4-branes and D2-branes in (4.6). 
The theory on D4-branes is simply the 4d gauge theory on M 4 , and describing its IR physics 
via its M-theory lift was one of the main points of [20]. The theory on D2-branes is a 2d 
theory with J\f = (2, 2) supersymmetry preserved by the system (4.6), see e.g. [9,42-44]. This 
2d theory couples to 4d gauge theory and, hence, describes a half-BPS surface operator as a 
combined 2d-4d system. 

This has to be compared with our earlier discussion in Section 3.2, where we saw that sur¬ 
face operators constructed from codimension -2 defects naturally lead to singularities of gauge 
fields in the 4d gauge theory, while now we see that surface operators built from codimension-4 
defects naturally lead to a description via combined 2d-4d system. Furthermore, the number 
N of D4-branes that determines the rank of the gauge group in four dimensions is the rank of 
the flavor symmetry group from the viewpoint of 2d theory on the D2-branes. In particular, 
in the basic case of = 2 each D2-branes carries a 1/(1) linear sigma-model with N = 2 
charged flavors, whose Higgs branch is simply the Kahler quotient C‘^//U{1) = CP^. 

This implies that codimension-4 defects give rise to a 2d-4d coupled system, in which 
gauge theory in the bulk is coupled to the CP^ 2d sigma-model on D C M 4 , which is IR- 
equivalent to the corresponding 2d gauged linear sigma model. Moreover, this also shows why 
the deformation associated to the separation along direction in (4.5) does not affect the 
corresponding twisted superpotential. And here the identification of unbroken supersymmetry 
and the precise type of the F-terms in 2d becomes crucial. 

Namely, from the viewpoint of the D2-branes in (4.6), the separation along the direc¬ 
tion is the gauge coupling constant of the 2d gauged linear sigma-model [9,42-44], 




52d — 




(4.12) 


D2 


On the other hand, it is a standard fact about 2d Af = (2, 2) supersymmetry algebra that 
twisted superpotential is independent on the 2d gauge coupling constant [45]. 

The reader may observe that the number of variables Ui parametrizing the positions of 
the created M2-branes exceeds the number of parameters introduced via (3.3) for surfaces 
of genus g > 1. At the moment it does not seem to be known how exactly one may describe 
the system with M5- and Mh^branes at an intermediate energy scale in terms of a four¬ 
dimensional quantum field theory. It seems quite possible that the resulting description will 


involve coupling one gauge field 
D C M4, in general. 


(r) 


to more than one copy of the CP^ 2 d sigma-model on 


4.3 Twisted superpotentials as generating functions 

As we have seen in the previous subsection, regardless how different the theories with two 
types of surface operators may be in the UV, their effective descriptions in the IR have a 
relatively simple and uniform description. More specifically, the theories we are considering 
in this paper are essentially determined in the IR by their twisted superpotentials. Hence we 
focus on them. 


- 24 - 




The twisted superpotentials in the presence of codimension-2 and codimension-4 surface 
operators will be denoted by and , respectively. The twisted superpotential = 
W^®(a, X, t) depends besides a and r on coordinates x for BunG(C), and u, r) 

on the positions of the points on C where the codimension-2 defects are located. 

From both W”^® and we can find the corresponding Gukov-Witten parameters 

t'^®(a, x,r) and r) via (4.10). If the two surface operators are equivalent in the 

deep IR there must in particular exist an analytic, locally invertible change of variables 
u = u^{x-, a, t) relating the Gukov-Witten parameters t and t' as 


X, t) = a, r), r). 


(4.13) 


It follows that the twisted superpotentials W^® and may differ only by a function 

independent of o. 

One may furthermore note that the variables ut are dynamical at intermediate scales, or 
with non-vanishing Omega-deformation. The system obtained by separating the M5'-branes 
by some finite distance from the M5-branes will be characterized by a superpotential 
W' depending both on x and u, in general. We had argued above that this superpotential 
does not depend on the separation Ax^. Flowing deep into the IR region one expects to 
reach an effective description in which extremization of the superpotential determines u as 
function of x and the remaining parameters, u = m*(x, a, r). The result should coincide with 
VV”^®(a, X, r), which is possible if the resulting superpotential W' differs from W’^^(a, u, r) by 
addition of a function W" {u,x,t) that is a-independent 

>V'(a, x,rt, r) = >V“^(a, u, r)-|->V^'(u, x, r); (4-14) 


the additional piece W” {u, x,t) may be attributed to the process creating the M2-branes 
from M5hbranes. Extremization of W' implies that 




= -j^>V"(u,X,t) 

u=u*[x,a,T) O'U'f 


w=w*(ir,a,r) 


(4.15) 


and VV’'(a, X, tt, r)|^^^ should coincide with VV’'^®(a, x, r). 

We are now going to argue that W'^®, and W represent generating functions for 
changes of variables relating three different sets of Darboux-coordinates for the same moduli 
space A42d locally parametrized by the variables a and x (see, for example, [46], Section 2.1, 
for the definition of generating functions and a discussion of their role in the Lagrangian 
formalism). 

Considering W^® first, one may define other local coordinates for A42d as 

Pr = -^>V“®(a,x,T). (4.16) 


Both (x,p) and (a,t), with t defined via (4.10), will generically define local coordinates for 
A42d- Having a Poisson-structure on A42d that makes (x,p) into Darboux-coordinates it 
follows from (4.10) and (4.16) that (a,t) will also be Darboux-coordinates for A42d- 
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If X and u are related by a locally invertible change of variables u = a, r) it follows 
from (4.15) that u together with the coordinates v defined by 

Vr = , (4.17) 

OUr 

will represent yet another set of Darboux coordinates for Ai 2 d- In this way one may identify 
and W as the generating functions for changes of Darboux variables (a, t) •(-)• («, v) and 
{u,v) •(-)• {x,p) for Alvac, respectively. 

There are various ways to compute the twisted superpotential W. One (though not the 
only one!) way is to compute the asymptotic expansion of the Nekrasov partition function [26] 
in the limit ei ^2 0. It takes the form 


F W 

log + (4.18) 

6162 61 

Here, F is the Seiberg-Witten prepotential that does not depend on the surface operator and 
defines the corresponding IR 4d theory in the bulk. The next term in the expansion, W, is 
what determines the IR theory with the surface operator.^ 

In what follows we will use the relations of the instanton partition functions to conformal 
blocks to determine W^®(a, x, r) and VV’'^^(a, u, r) via (4.18). Both functions will be identified 
as generating functions for changes of Darboux-variables {x,p) -H- {a,t) and {u,v) -H- {a,t) 
for the Hitchin moduli space A^h(C'), respectively. Among other things, this will imply that 
VV”^®(a, x, t) and VV”^^(a, m; u) indeed satisfy a relation of the form 


W^®(a,x,T) = W'^^(o, rt*(x, a, r), r) -I- W®°'^(u*(x, a, r), X, r). 


(4.19) 


In view of the discussion above one may view this result as nontrivial support for the conjec¬ 
tured IR duality relation between the theories with the surface operators of co-dimensions 2 
and 4, if we set W" = 


4.4 Relation to conformal field theory 

We had previously observed that the twisted superpotentials VV’^®(a, x,r) that may be cal¬ 
culated from the instanton partition functions Z^®(a, x, r; ei, € 2 ) via (3.15) represent changes 
of Darboux variables for the Hitchin integrable system. We will now discuss analogous re¬ 
sults for >V^^(a, u, r). To this aim we begin by describing the expected relations between the 
instanton partition functions Z^^{a,x,T; € 1 , 62 ) and Liouville conformal blocks. 

Conformal blocks for the Virasoro algebra with central charge Cb = 1 -\- 6(6 -|- 6~^)^ 
may be defined in close analogy to the Kac-Moody conformal blocks discussed above. Our 
discussion shall therefore be brief. Given a Riemann surface C with n punctures, we associate 
representations Va^ generated from highest weight vectors Va^ to the punctures Zr, r = 1,... ,1. 
The Lie algebra Vect(C'\{ 2 :i,..., zi}) of meromorphic vector fields on C with poles only at Zr, 

a system without surface operators one has VV = 0. 
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r = 1,is naturally embedded into the direct sum of I copies of the Virasoro algebra with 
the central elements identified (using the expansion of the vector fields near the punctures). 
Conformal blocks ip are then defined as linear functionals on Va^ that are invariant 

under the action of Vect(C' \ {zi, ..., zi}). This invariance condition represents the conformal 
Ward identities. Chiral partition functions Zjr(ip,C',b) are defined as the evaluation of ip on 
the product of highest weight vectors 0^=1 in the physics literature often denoted as 

Z^{ip, C-b) = { ^ . (4.20) 

In general, the space of conformal blocks is infinite-dimensional. However, it can be 
decomposed into a direct sum (or direct integral, depending on the situation) of finite¬ 
dimensional spaces (in some cases, such as that of the Liouville model, one-dimensional 
spaces, so that we obtain a basis) using the gluing construction reconstructing C from its 
pants decompositions specified by the data a = (C, T) introduced in Section 2.1. Its elements 
are labeled by representation parameters j5e assigned to the cut curves je G C. We denote 
the resulting chiral partition functions by Z^(l3,T;b). 

We shall also discuss the situation of d additional degenerate representations V_i/ 2 fe 
(sometimes called $ 1^2 primary fields) associated to points S = {ui,... ,U(i} C C that are 
distinct and different from the punctures zi,... ,zi. The corresponding chiral partition func¬ 
tions then satisfy d second order differential equations resulting from the existence of degree 
2 null vectors in V_i/ 2 b- A basis for the space of solutions can be obtained by starting from 
a pants decomposition a of C. Each pair of pants Cq 3 obtained by cutting along C contains 
a subset Sv of S. Choosing a pants decomposition of Cg 3 \ Sv one obtains a refined pants 
decomposition a that can be used to define chiral partition functions Z^ ^{(5, u, q] b) as before. 
The additional set of labels w entering the definition of Z^ ^ is constrained by the fusion rules 
for existence of conformal blocks with degenerate representations inserted, and may therefore 
be represented by elements of Z^. 

The precise definition of the instanton partition functions Z^'^ = Z^'^ in the presence 
of d codimension 4 surface operators depends on the choice of a refined pants decomposition 
a, decorated with certain additional discrete data collectively denoted vu, see [44]. In [9] it 
was conjectured that the instanton partition functions Z^'^ coincide with Liouville conformal 
blocks with d additional degenerate helds inserted, 


2^S^mia-,u,T-,ei,e2) = Z^^^{f3,u,T]b), 
given that the parameters are related as 


o _ Q . ■ 

Pe — t: + 1 . - 5 

2 y/eie2 


b^ = 

£2 


(4.21) 


(4.22) 


Further evidence for (4.21) and some of its generalizations were discussed in [6,44,47,48]. 


Now we are ready to bring together the results of the previous sections to demonstrate 
the IR duality of two 4d gauge theories with surface operators and to link it to the separation 
of variables in CFT and Hitchin system. 
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4.5 Relation to the Hitchin system and to the separation of variables 

It is shown in the Appendix D that (4.21) implies that 


logZ^^{a,u,T-ei,€2) ~- —T{a,T)-—W^^{a,u,T), 

eie2 €i 


(4.23) 


as already proposed in [9]. The function VV’^^(a,rt, r) is given as 

h 


/ Uff. 

V. 

K=L 


(4.24) 


We are now going to explain that there exist other sets of natural Darboux-coordinates {u, v) 
for Hitchin moduli space allowing us to identify the function W’'^^(a, u,r) defined in (4.24) 
as the generating function for the change of variables (a, t) •(-)• {u,v). 

Recall from Section 2.3 that the spectral cover construction allows us to describe AIh(C') 
as the space of pairs (S,L). The line bundle L may be characterized by a divisor of zeros 
of a particular section of L representing a suitably normalized eigenvector of the Higgs field 
if G H^{C, End(T) ( 8 > Kc) that we describe presently. Even though this divisor is not unique, 
it’s projection onto C is uniquely determined by the data of the rank two bundle B with a 
fixed determinant^ and the Higgs field ip. 

Locally on (7, we can trivialize the bundle B and choose a local coordinate z. Then we 
can write (p as 


(p = 


faiz) 

Vc(2) 


biz) \ 

-a(z)J 


dz. 


We have the following explicit formula for the eigenvectors of (p 




(aiv) ± viy)\ 

V c(z) 


v'^iy) 


^tT{(p‘^{y)). 


Note that for the matrix element ciz)dz to be well-defined globally on C and independent of 
any choices, we need to represent B as an extension of two line bundles, see Appendix C.2 
for more details. 

If c(z) 7 ^ 0, then T 7 ^ 0 for either branch of the square root. If c(z) = 0, then one of them 
vanishes. Now recall that the line bundle L on the double cover S of C is defined precisely 
as the line bundle spanned by eigenvectors of p) (at a generic point p oi C, p has two distinct 
eigenvalues, which correspond to the two points, p' and p", of S that project onto p, and the 
fibers of L over p' and p” are the corresponding eigenvectors). Therefore, if we denote by 
V the divisor of zeros of c{z)dz on C, T gives rise to a non-zero section of L outside of the 
preimage of P in S. 

®As explained in Appendix C.2, a natural possibility is to consider rank two bundles B whose determinant 
is a fixed line bundle of degree 2g — 2 + n. The moduli space of such bundles is isomorphic to the moduli space 
of SL 2 -bundles on C. 


- 28 - 




Generically, T) is multiplicity-free and hence may be represented by a collection u = 
{ui ,..., Ud) of d := deg(P) distinct points. The number number d depends on the degrees 
of the line bundles used to represent B as an extension, in general. It may be larger than 
3g — 3 + n, the dimension of Bunc*. However, fixing the determinant of B defines a collection 
of constraints allowing us to determine Uk, k = h + 1,d in terms of the coordinates Ui, 
i = l,...,Uh. 

There are two distinct points, u' and u'l, in S over each Ui £ C. Then for each i = 1,..., h, 
our section has a non-zero value at one of the points, u'^ or u'l, and vanishes at another 
point. Thus, the divisor of this section on S is the sum of particular preimage of the points 
Ui,i = 1,... ,h, in E, one for each i. While there is a finite ambiguity remaining for this 
divisor,® the unordered collection u = {ui ,..., Uh) of points of C is well-defined (generically). 
And then for each Ui we choose the eigenvalue £ T^C, for which our section provides a 
non-zero eigenvector. It is known that the collection {u,v) = ((tti,ui),..., {uh,Vh)) can be 
used to get to a system of Darboux coordinates for Aiu{C) [49,50], see also [51] for related 
results. 

It was observed in [50] that the definition of the variables {u, v) outlined above can be seen 
as a generalization of the method called separation of variables in the literature on integrable 
models [16]. A familiar example is the so-called Gaudin-model which can be identified with the 
Hitchin integrable system associated to surfaces C of genus zero with n regular singularities 
at distinct points zi,..., Zn- The Higgs field can then be represented explicitly as 





^ ^ Ar — 0, 

r=l 


where 




Mr A 
Ur- MV’ 


and the separated variables are obtained as the zeros of the lower left entry A {y)dy of (p: 

UkZfiy-Uk) 


A (y) = u 

n—1 




Vk = 


E 

r=l 


A® 


Uk Zr 


(4.25a) 

(4.25b) 


One may think of the separation of variables as a useful intermediate step in the construc¬ 
tion of the mapping from the original formulation of an integrable model to the description 
as the Hitchin fibration in terms of action-angle coordinates (a,t). The remaining step from 
the separated variables (n, v) to the action-angle variables is then provided by the Abel map. 
The function VV”^^(a, tt, r) is nothing but the generating function for the change of Darboux 
coordinates between {u,v) and (a, t). A few more details can be found in Appendix C.4. 

®More precisely, we have choices of the preimages u'i or u” for each i, which agrees with the number 

of points in a generic Hitchin fiber corresponding to a fixed SL 2 bundle. 
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4.6 IR duality of surface operators from the defects of codimension 2 and 4 

In this section we combine the ingredients of the brane analysis in Section 4.1 with our results 
on the twisted superpotentials to show that the 4d gauge theories with the surface operators 
constructed from codimension-2 and codimension-4 defects are equivalent in the IR. 

Indeed, their vacuum structures are controlled by the twisted superpotentials W^®(a, x, r) 
and tt, r), and we have found that they are related by a change of variables (that is, 

a redefinition of fields). 

Furthermore, when combined, the above arguments - including the brane creation upon 
the change of separation in the direction - show that two types of surface operators 
constructed from codimension-2 and codimension-4 defects preserve the same supersymmetry 
subalgebra and have the same twisted chiral rings.^ This is sufficient to establish their 
equivalence for the purposes of instanton counting. In order to demonstrate the IR equivalence 
of the full physical theories, we need to show the isomorphism between their chiral rings (and 
not just the twisted chiral rings). In general, this is not guaranteed by the arguments we 
have used, but the good news is that for simple types of surface operators, including the 
ones considered here, the chiral rings are in fact trivial® and, therefore, we do obtain the 
equivalence of the two full physical theories. 

As we already mentioned in the Introduction, this equivalence, or duality, between the IR 
physics of 4d AA = 2 gauge theories with two types of surface operators is conceptually similar 
to the Seiberg duality of 4d A/" = 1 gauge theories [1]. In fact, it would not be surprising if 
there were a more direct connection between the two phenomena since they both enjoy the 
same amount of supersymmetry and in its brane realization, Seiberg’s duality involves the 
same kind of “moves” as the ones described in the previous section. 

4.7 Turning on the Omega-deformation 

The relation between W”^®(a, x,r) and W'^^(a,«, r) has a rather nontrivial generalization in 
the case of non-vanishing Omega-deformation that we will describe in this subsection. The 
fact that in 2d this a variant to the separation of variables continues to hold for non-zero 
values of ei and €2 suggests that the two 4d AA = 2 gauge theories remain IR equivalent even 
after Omega-deformation. The possibility of such an equivalence certainly deserves further 
study. 

When we quantize the Hitchin system, the separation of variables may also be quantized. 
In the genus zero case, in which the quantum Hitchin system is known as the Gaudin model, 
this was first shown by E. Sklyanin [16]. Note that the quantization of the classical Hitchin 
system corresponds, from the 4d point of view, to “turning on” one of the parameters of the 

^Twisted chiral rings are Jacobi rings of the twisted chiral superpotential VV which has been our main 
subject of discussion in earlier sections. 

®In general, 2dAf= (2, 2) theories may have non-trivial chiral and twisted chiral rings, see for example [52]. 
However, if we start with a 2d theory without superpotential, then, as long as chiral superfields are all massive 
in the IR, integrating them out leads to a theory of twisted chiral superfields with a twisted superpotential, 
and so the chiral ring is indeed trivial. 
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Omega-deformation which is the case studied in [53]. It has been explained in Section 6 of [17] 
that one may interpret the separation of variables in the Gaudin model, as well as more general 
quantum Hitchin systems, as the equivalence of two constructions of the geometric Langlands 
correspondence (Drinfeld’s “first construction” and the Beilinson-Drinfeld construction). 

Feigin, Frenkel, and Stoyanovsky have shown (see [18]) that in genus zero the separation 
of variables of the quantum Hitchin system maybe further deformed when we “turn on” 
both parameters of the Omega deformation. This result was subsequently generalized to 
get relations between non-chiral correlation functions of the WZW-model and the Liouville 
theory in genus 0 [19], and in higher genus [54]. It has furthermore been extended in [27] 
to larger classes of conformal blocks. From the 4d point of view, this relation amounts to a 
rather non-trivial relation via an integral transform (a kind of “Fourier transform”) between 
the instanton partition functions of the Omega-deformed 4d theories with surface operators 
corresponding to the defects of codimensions 2 and 4. 

The resulting relation has its roots in the quantum Drinfeld-Sokolov reduction. We 
recall [55,56] that locally it amounts to imposing the constraint J~{z) = 1 on one of the 
nilpotent currents of the affine Kac-Moody algebra s[ 2 . The resulting chiral (or vertex) 
algebra is the Virasoro algebra. Furthermore, if the level of s [2 is 

^ = - 2 - 

then the central charge of the Virasoro algebra is 

c = 1 -|- 6(6 -|- 6“^)^. 

Globally, on a Riemann surface C, the constraint takes the form J~{z)dz = uj, where cu is 
a one-form, if we consider the trivial S'L 2 -bundle, or a section of a line bundle if we consider a 
non-trivial 5'L2-bundle that is an extension of two line sub-bundles (the representation as an 
extension is necessary in order to specify globally and unambiguously the current J~{z)dz). 
Generically, oj has simple zeros, which leads to the insertion at those points of the degenerate 
fields R_i/ 2 b of the Virasoro algebra in the conformal blocks. 

It is important to remember that classically the separated variables Ui are the zeros of a 
particular component of the Higgs field if. But the Higgs fields correspond to the cotangent 
directions on parametrized by the p-variables. After quantization, these variables 

are realized as the derivatives of the coordinates along the moduli of 5L2-bundles (the x- 
variables), so we cannot directly impose this vanishing condition. Therefore, in order to 
define the separated variables u in the quantum case, we must first apply the Fourier transform 
making the p-variables into functions rather than derivatives (this is already needed at the 
level of the quantum Hitchin system, see [17]). Since the Fourier transform is an integral 
transform, our formulas below involve integration. Indeed, the separation of variables linking 
the chiral partition functions in the WZW-model and the Liouville model is an integral 
transform. 
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In Appendix E it is shown that the relations described above can be used to derive the 
following explicit integral transformation, 


= Nj / dui . . . dUn-3 {x,u) Z^{u, z) 


(4.26) 


where Z^"^ and ZZ are obtained from Z^"^ and Z^ by taking the limit Zn ^ oo, and the 
kernel ICsoy{x,u) is defined as 




' n—1 

E 


r=l 


, mZliZr-U,) 


n—S n—1 


k<l 


r=l 


' u:^^izr-z.y 

nlZliZr - U,) 


Nj is an (x, 2 ;)-independent normalization factor that will not be needed in the following. 
Note that the x-dependence it entirely in the hrst factor on the right hand side of (E.12). 
Using (3.15), (4.23) and (E.14) it is easy to see that the relation (4.19) follows from (4.26). 
Eormula (4.26) is the relation (1.1) discussed in the Introduction made explicit. 

Thus, we see that the separation of variables in the most general case (with both parame¬ 
ters of the Omega deformation being non-zero), viewed as a relation between the chiral chiral 
partition functions in the WZW-model and the Liouville model, provides the most satisfying 
conceptual explanation of the IR duality of the 4d gauge theories with surface operators of 
two kinds discussed in this paper. 


APPENDICES 

A. Surface operators and Nahm poles 

Complex (co)adjoint orbits are ubiquitous in the study of both half-BPS surface operators 
and boundary conditions. This happens for a good reason, and here we present a simple 
intuitive explanation of this fact. In short, it’s due to the fact that both half-BPS surface 
operators and boundary conditions are labeled by solutions to Nahm equations. Then, the 
celebrated work of Kronheimer [57] relates the latter to complex coadjoint orbits. 

Suppose that in our setup (3.1) we take C = x M and M 4 = D x D' = where 
D' = is the “cigar.” In other words, D' is a circle fibration over the half-line, M+ = {y > 0}, 
with a singular fiber at y = 0 so that asymptotically (for y —)■ -|-cx)) D looks like a cylinder, 
see Figure 3. Then, the six-dimensional (2, 0) theory on M 4 x C with a codimension-2 defect 
on T) X C can be reduced to five-dimensional super-Yang-Mills theory in two different ways. 
First, if we reduce on a circle S^, we obtain a 5d super-Yang-Mills on M 4 x M = with a 
surface operator supported on H x M = M^. If we denote by r = the radial coordinate 
in the plane transverse to the surface operator, then the supersymmetry equations take the 
form of Nahm’s equations: 

da r, 1 db , , dc . .n 

- = [ 6 ,c], ^ = [c,«l, ^ = 1 “.'']. (A.l) 
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Figure 3: The six-dimensional (2,0) theory with a codimension-2 defect at the tip of the cigar reduces 
to 5d super-Yang-Mills theory with a non-trivial boundary condition. 

where we used the following ansatz for the gauge field and for the Higgs field: 


A = a{r)d6 . 


uT 

(j) = b{r) - \-c{r)d9 . 

r 


On the other hand, if we first reduce on the circle fiber Sp of the cigar geometry D', we 
obtain a 5d super-Yang-Mills on M_|_ x D x C with a non-trivial boundary conditions at 
y = 0 determined by the codimension-2 defect of the six-dimensional theory. Note, these 
boundary conditions are also associated with solutions to Nahm’s equations (A.l) for the 
Higgs field (p = (a, 6, c). Further dimensional reductions of these two systems yield many half- 
BPS boundary conditions and surface operators in lower-dimensional theories, all labeled by 
solutions to Nahm’s equations. 

Among other things, this duality implies that 
similar physical and mathematical structures can 
be found on surface operators as well as in the 
study of boundaries and interfaces. A prominent 
example of such structure is the algebra of pa¬ 
rameter walls and interfaces, i.e. Janus-like soli- 
tons realized by monodromies in the space of pa¬ 
rameters. (In the case of surface operators, such 
monodromy interfaces are simply line operators, 
which in general form non-commutative algebra 
if they can’t move off the surface operator, as 
illustrated in Figure 4.) 

This description of walls, lines and interfaces 
as monodromies in the parameter space provides 
a simple and intuitive way of understanding their 

non-commutative structure and commutation relations; it is captured by the fundamental 
group of the parameter space [7]: 



Figure 4: In the presence of surface operator 
and/or Omega-background line operators do 
not commute. 


TTi ({parameters}) 


(A.2) 
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For instance, in the case of C = one finds tti {{Tc/S nY^^), which is precisely the braid 
group (in the case, of type Aj\f-i). It is generated by parameter walls / interfaces Li that 
obey the standard braid group relations: 

Li -k Li-i-i -k Li = Lj+i k Lik Li-i-i (A-3) 

From 2d and 3d perspectives, these systems are often described by sigma-models based on 
flag target manifolds (or their cotangent bundles) where the lines/walls Li are represented 
by twist functors; see [8,41] for further details and many concrete examples of braid group 
actions on boundary conditions. The case of the parameter space (2.5) is qualitatively similar. 

B. Twisting of Kac-Moody conformal blocks 

This appendix collects some relevant mathematical background concerning the dependence 
of Kac-Moody conformal blocks on the choice of a holomorphic bundle on C. 

B.l Twisted conformal blocks 

A generalization of the defining invariance condition allows us to define a generalized notion 
of conformal blocks depending on the choice of a holomorphic G-bundle B on C. One may 
modify the defining invariance condition (3.7) by replacing the elements of the Lie algebra 
flout by a section of 

flout := r(C,fle), flg:=exg. (B.l) 

(j 

Describing B in terms of a cover {77*; z G X} of C allows us to describe B in terms of the 
G-valued transition functions hij{z) defined on the intersections 77*j = 77* 0 77^. The sections 
of flout are represented by families of g-valued functions 77 * in 77*, with ry* and rjj related on 
the intersections 77*^ by conjugation with hij{z). In this way one defines B-twisted conformal 
blocks ifQ depending on the choice of a G-bundle B. 

More concrete ways of describing the twisting of conformal blocks are obtained by choos¬ 
ing convenient covers {77*; z G X}. One convenient choice is the following: Let us choose discs 
Dfc around the points Zk, k = 1,... ,n such that 77out ■= C \ {zi, ..., z**} and 77in = Ufc=i 
form a cover of C. It is known that for G = SL{2) G-bundles B can always be trivialized 
in 77out and 77in. An arbitrary G-bundle B can then be represented by the G-valued tran¬ 
sition functions h^^tk) defined in the annular regions := 77out O D/* modulo changes of 
trivialization in 77in and in 77out, respectively. 

Introducing the dependence on the choice of B in the way described above makes it easy 
to see that infinitesimal variations 6 oi B can be represented by elements of g ® C((7*)). 
Choosing a lift Xs to the diagonal central extension of g ® C((7*)) allows us to define a 
(projective) action of TBunclg on CBg(77.i,..., TZn)- This means that a differential operator 
5 representing an element TBunclg can be represented on the conformal blocks in terms of 
the action of qs on (S>r=i'^r, schematically 

<JA(e[n]) = , e[**] := ei ® • • • (g) e„ . (B.2) 
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This action describes the response of a conformai biock ipQ with respect to an inhnitesimai 
variation of B. 


B.2 Genus zero case 

In the case of genus 0 it suffices to choose the transition functions hk{tk) in the annuiar regions 
Afc around the points to be the constant niipotent matrices = (g ) • The coiiection 

of parameters x = (xi,..., Xn) can be used to represent the dependence on the choice of B in 
this case. The action of TBunclg on spaces of conformai biocks defined via (B.2) may then 
be represented more expiicitiy in terms of the differentiai operators defined as 

Jr" = dxr, Jr = Xrdxr “ Jr, = -xldx, + 2jr®r • (B.3) 


The Casimir operator is represented as muitipiication by jr{jr + !)• 

The parametrization in terms of n variabies x = (xi,...,x„) is of course redundant. 
The conformai Ward-identities (3.7) inciude the invariance under giobai s[ 2 -transformations, 
aiiowing us to eiiminate three out of the n variabies xi,..., in the usuai way. 

The operators Hr appearing in the Knizhnik-Zamoiodchikov equations (3.11) are then 
given by the formulae 


Hr 



(B.4) 


where the differential operator J^s is defined as 


Jrs := W Jr“ J/ := JrJs + JJ + JJ Js) 


(B.5) 


The operators Hr commute, and may therefore be used as Hamiltonians for generalizations 
of the Gaudin models associated to more general representations of SL{2,C). 


B.3 Higher genus cases 

Instead of the covers considered in Subsection B.l above one may use alternatively use covers 
defined using the gluing construction. One thereby gets a cover {Ui; i €l} with intersections 
represented by annuli Ag between pairs of pants or connecting two legs of the same pair of 
pants. Choosing constant diagonal transition functions ^-i) in the annuli Ag gives us a 
collection of local coordinates Xg, e = 1,..., 3^ — 3 + n for Bunc, G = SL{2). The resulting 
parameters x for Bunc are easily identified with the parameters x introduced in the gluing 
construction of conformal blocks via (3.10) provided we choose /C(r, x) to be 

In order to have a globally well-defined current J~ on C one needs to represent B as an 
extension. Taking 

(B.6) 

appears to be particularly natural. This allows us to represent J~ as a section of C®Kc- As 
explained in Appendix C it is natural in our case to consider fixed line bundles C of degree 
d!. Let us represent C as 0{V'), with divisor V being represented by the points yi,.. .yd'- 
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The bundle 13 may be described by using a cover {Ui \ ^ G X} for C containing small discs 
around k = 1,... with transition functions 


h'k 



(B.7) 


on the annuli = D^\{yfc}, where is a coordinate on vanishing at y^. Sections of B may 
alternatively be represented locally by functions that are regular outside oi {yk,k = 1,... .,d'} 
and may have poles with residue in a fixed line at y^, k = 1,... ,d'. Using the transition 
functions (B.7) determines the lines ik in terms of the parameters Xk. Modifications of B that 
increase the degree d' of X are called Hecke modifications. 

Using covers defined with the help of the gluing construction it appears to be natural to 
take d' = 2g — 2. In this case one may assume that there is exactly one yk contained in each 
pair of pants. Kac-Moody conformal blocks associated to each pairs of pants appearing in the 
pants decomposition of a closed Riemann surface can then be defined using conformal blocks 
on with one insertion being the degenerate representation of the Kac-Moody algebra 
77fc/2 representing the Hecke modifications within conformal held theory [27]. If the Riemann 
surface has punctures, one may use conformal blocks on Co ,3 without extra insertion of 'R-k /2 
for the pairs of pants containing the punctures. 

It is worth remarking that d! = 2g — 2 \s exactly the case where the current J ~, being a 
section of Kq 0 X, has 4^ — 4 zeros Ui, as required by the identihcation of the points Ui with 
the end-points of the M2-branes created from the M5'-branes. 


C. Holomorphic pictures for the Hitchin moduli spaces 

The Hitchin space A4 h(C') was introduced in the main text as the space of pairs {B,ip). 
Interpreting the Higgs helds ip G H^{C, End(X) 0 Kc) as representatives of cotangent vectors 
to Bunc, one may identify A4 h(C') with T^Bunc, the cotangent bundle of the moduli space 
of holomorphic G-bundles on C. This description equips A4 h(C') with natural complex and 
symplectic structures, leading to the definition of local sets of Darboux coordinates {x,p) 
parametrizing the choices of G-bundles via coordinates x, and the choices of Higgs fields (p in 
terms of holomorphic coordinates p. 

In order to exhibit the relation with conformal field theory we will find it, following [27,58], 
useful to consider a family of other models for AIh(G). We will consider moduli spaces M^{C) 
of pairs (H, V^) consisting of holomorphic bundles B with holomorphic e-connections V^. An 
e-connection is locally represented by a differential operator = {edy+A{y))dy transforming 
as Vg = g~^ • Vg - 5 under gauge-transformations. Consideration of M^{C) will represent a 
useful intermediate step which helps clarifying the link between conformal field theory and 
the Hitchin system. Noting that any two e-connections and differ by an element of 
77'^(C', End(X) ® Kq) one sees that M^{C) can be regarded as a twisted cotangent bundle 
TgBunG. Picking a reference connection V^q, one may represent a generic connection as 
V'=V',o + </^. 
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To avoid confusion let us stress that the resulting isomorphism 7Wg(C') ~ T^Bunc is not 
canonical, being dependent on the choice of V'^q. Instead we could use the known results of 
Hitchin, Donaldson, Corlette and Simpson [59-63] relating pairs {B, ip) to flat connections on 
C to identify the moduli spaces M.^{C) and 7W^(C'). The description of as twisted 

cotangent bundle yields natural complex and symplectic structures which are inequivalent for 
different values of e. This can be used to describe the hyperkahler structure on A^h(C'), with 
e being the hyperkahler parameter [64], 

However, in order to discuss the relation with conformal field theory we find it useful to 
adopt a different point of view. The definition of conformal blocks depends on the choice of 
a G-bundle B, which may be parametrized by variables x in a way that does not depend on 
ei and e 2 - The gluing construction yields natural choices for the reference connection V^q, 
e.g. the trivial one. All dependence on the parameter e is thereby shifted into the relations 
between different charts 14^ on AI^(C') parametrized in terms of local coordinates in 

a way that does not explicitly depend on e. 

One may formally identify G End(T) ® Kq) as an e-connection for e = 0. We 

therefore expect that the Darboux coordinates (xe,pe) turn into the Darboux coordinates 
{x,p) discussed in the main text when e —)■ 0. This will be further discussed below, after 
having discussed possible choices of Darboux coordinates more concretely. 

C.l Three models for Hitchin moduli space 

There are three models for of interest for us: 

(A) As space of representations of the fundamental group 

Hom(7ri(C), SL(2, C))/SL(2, C). (C.l) 

(B) As space of bundles with connections (T, V^), 

V; = (td, + .4(rf) dy , A{y) = ( ) , (C.2) 

Having n punctures zi,... ,Zn means that A{y) is allowed to have regular singularities 
at y = Zr of the form 

^{y) = - - -h regular. (C.3) 

y- Zr 

(B’) As space of opers e^dy -|- t{y), where t{y) has n regular singularities at y = Zr, 

t{y) = -—-regular, (C.4) 

[y - Zry y- Zr 

and d apparent singularities at y = Uk, 

t{y) = ^ + regular , (C.5) 

4(y - UkY y-Uk 
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Having an apparent singularity ai y = Uk means that the monodromy around is 
trivial in PSL(2,C). This is known [17, Section 3.9] to be equivalent to the fact that 
the residues Hr, r = 1,... ,n are constrained by the linear equations 

vl + tk ,2 = 0, k = l,...,l, t{y) = '^tk,i{y-Uky~‘^. (C.6a) 

1=0 

If g = 0, the parameters Hg, s = 1,..., n are furthermore constrained by 

n 

ZriZrHr + {a + l)6r) = 0 , o = —1,0,1, (C.6b) 

r=l 

ensuring regularity of t{y) at infinity. 

Models (B) and (B’) are related by singular gauge transformations which transform A{y) to 
the form 

(C.7) 

In order to describe the relation between (B) and (B’) more concretely let us, without loss of 
generality, assume that elements of Bunc are represented as extensions 

0 ^ ^ B ^ C" ^0. (C.8) 


Describing the bundles B by means of a covering Ui of C and transition functions Bij between 
patches Ui and Up one may assume that all Tjj are upper triangular. 


Bij 


0 \ A TA 
V 0 c'j Vo 1; • 


(C.9) 


This implies that the lower left matrix element A~{y) of the e-connection edy A{y) is a 
section of the line bundle {C)~^ ® C" ® Kc, with Kq being the canonical line bundle. The 
gauge transformation which transforms A{y) to the form (C.7) will be singular at the zeros 
Uk of A~{y), leading to the appearance of the apparent singularities Uk in (C.4). 


C.2 Complex-structure dependent Darboux coordinates 

Let us briefly discuss possible ways to introduce Darboux coordinates {x,p) for and 

how the passage from e-connections to opers dehnes a change of Darboux coordinates from 
(x,p) to {u, v). 


Genus zero 

In the cases of genus 5 = 0 
Ar — 


we may parametrize the matrices Ar in (C.2) as 

/ A^ T+ \ _ A —Ar- ^ ^ A 

-HO j ^ Vo 1 y \Pr -Ir) Vo 1 y ’ 


(C.IO) 
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assuming that {xr,Pr) are a set of Darboux coordinates with {pr,Xs} = Sr,s- Let Vn be the 
phase space whose algebra of functions is generated by functions of {xr,Pr), r = 1,... ,re. 
The space A^flat(C'o,n) can be described as the symplectic reduction of Vn w.r.t. the global 
sl 2 -constraints 

n 

Y,^r=0, (C.ll) 

r=l 

for a = —, 0 ,+, or, more conveniently, as the symplectic reduction of Vn-i w.r.t. the con¬ 
straints (C.ll) for a = —, 0 combined with sending Zn —)• oo. We will use the latter description. 

The change of {x,p) •(-)• {u,v) induced by the relation between models (B) and (B’) is 
explicitly described by the formulas (note that the same formulas (4.25) appear in the limit 
e ^ 0): 


A (y) 


Kzliy-zr) ’ 


Vk := A^{uk ), 


n—1 


^“(!/) = E 


y - Zr' 


(C.12a) 

(C.12b) 


The resulting change of variables (x,p) {u,v) is known to be a change of Darboux coor¬ 

dinates. It is in fact the classical version of the separation of variables transformation for 
the Schlesinger system [65]. In order to see this, let us consider in the model (B’) the case 
I = n — 3. In this case the equations (C. 6 ) determine the Hr as functions of the parameters 
{u, v), u = (ui,..., ui), V = {vi,... ,vi). The solutions Hr{u, v; z) to the constraints (C. 6 ) are 
the Hamiltonians of the Gamier system. The flows generated by the Hamiltonians Hr{u, v, z) 
preserve the monodromy of the oper e'^dy + t{y). 

In the model (B) one may consider the Schlesinger Hamiltonians defined as 


(C.13) 


It is well-known that the non-autonomous Hamiltonian flows generated by the Hr preserve 
the monodromy of the connection edy + A{y). The change of variables dehned via (C.12) 
relates the Hamiltonians Hr{x,p', z) to the Hamiltonians Hr{u, v, z) of the Gamier system. 


Higher genus 

Considering the cases of higher genus one may introduce Darboux coordinates associated to 
the model (B) as follows. To simplify the discussion slightly let us consider closed Riemann 
surfaces, n = 0. Representing the bundles B as extensions (C. 8 ), there are two places where 
the moduli may hide, in general: They may be hidden in the choice of the line bundles C', C'\ 
as well as in the extension classes S G H^{C~^), in terms of transition functions represented 
by the Tjj in (C.9). A particularly simple case is found by choosing C = O and L" = £ in 
(C. 8 ), with C being a fixed line bundle of degree 2g — 2. Fixing £ is equivalent to fixing the 
determinant of B. 
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The dimension of the space of extension classes is then dim(//^(£“^)) = g — l + deg(£) = 
3g — 3. The moduli of Bunc can therefore be parametrized by the choices of extension classes. 
Coordinates x = (xi,... on give coordinates for Bunc. 

Serre duality implies that the dual of is the space H^{C ® Kc)- Recall that 

the lower left matrix element A~ (y) of an e-connection edy + A{y) is a section of the line 
bundle C ® Kq- Finding coordinates for H^{C (8) Kc) that are dual to the coordinates x 
on with respect to the pairing provided by Serre duality will therefore give us 

coordinates p = {pi,... ,p 3 g- 3 ) that are canonically conjugate to the coordinates x on Bunc. 


C.3 Complex-structure independent Darboux coordinates 

Representing elements of in terms of the model (A) mentioned above allows one to 

introduce useful Darboux coordinates which do not depend on a choice of complex structure of 
C as opposed to the coordinates (rt, v) and (x,p) introduced before. A convenient description 
was given in [66] and references therein. 

Let us use the set-up from Section 2.1. A trivalent 
graph cj on C determines a pants decomposition de¬ 
fined by cutting along the simple closed curves 7 e which 
intersect the edge e of cr exactly once. For each (ori¬ 
ented) edge e we shall denote 7 e,s = 7e, 7e,t and ')e,u the 
simple closed curves which encircle the pairs of bound¬ 
ary components ( 7 e,l, 7 e, 2 ), (7e,2,7e,3) and (7e,l,7e,3), 
respectively, with labeling of boundary components in¬ 
troduced via Figure 5. Let Le,i := tr(p( 7 e,i)) for i G 
{s, t, u, 1, 2, 3,4}. One may represent Le,^, Le,t and Le,n 
in terms of Darboux coordinates a 
Poisson bracket 



and ke which have 


1 


Figure 5: Pants decomposition of 
four-holed sphere with a numbering of 
boundary components. 


{ Ue , k^i j- — 


(2vr) 


(C.14) 


The expressions are 


Le,s = 2 cosh(27rae/e) , (C.15a) 

Le,t{{Le,s)‘^ ~ ^) = 2(Le,2Le,3 + Ae,lLe,4) + Le,s(Ae,l£e,3 + £e,2Le,4) (C.15b) 

-b 2 COsh(27r/Ce/e) Ci2{Le,s)c34{Le,s) , 

Ae,u ((Le,s)^ ~ 4:) = 2(Le7Le,3 + £e,2Le,4) + Le,s(Re,2Le,3 + Re,lRe,4) (C.15c) 

-b 2cOsh(7r(2fce - ae)/e)-^Ci2(Le,s)c34(Le,s) , 

where Cij{Ls) is defined as 

Cig{Ls) = Ll + Lj + L] + LsLiLg - 4 . (C.16) 


Restricting these Darboux coordinates to the Teichmiiller component we recover the Fenchel- 
Nielsen length-twist coordinates well-known in hyperbolic geometry. 
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C.4 Limit e —)> 0: Recovering the Higgs pairs 

We now want to send e —>• 0. One may note that the equation {edy + A{y))'ip{y, x, z) can in 
the limit e be solved to leading order in e by an ansatz of the form 

= (C.17) 

where x(y; x, z) is an eigenvector of A{y) with eigenvalue v, 

My)x{y;x,z) = v{y)x{y,x,z). (C.i 8 ) 


The function v{y) representing the eigenvalue of A{y) must satisfy + t{y) = 0, where 


t{y) = 

Using t{y) we define the Seiberg-Witten curve as usual by 

S = {(v,u)\v'^ + t{u) = 0 } . 

Two linearly independent eigenvectors of A{y) are given by 


x±{y,x,z) = 


Ao{y) ±v 

A-iy) 


(C.19) 


(C.20) 


(C.21) 


One of x±{y] X, z) vanishes at the zeros of A-{y). It easily follows from these observations 
that the coordinates (x,p) and {u,v) for Ai^{C) turn into the coordinates for AIh(C') used 
in the main text when e —?• 0 . 

It follows from (C.17) that Oe and ke are in the limit €2 —)• 0 representable in terms of peri¬ 
ods of the canonical differential v on S. Given a canonical basis IB = {ai,..., a/i; , a))} 

for = 77i(S,Z)/77i(C',Z) one may dehne the corresponding periods as 


CLi — 




(C.22) 


For given pants decomposition a one may find a basis Bq- with the following property: For 
each edge e of cr there exists an index ie G {1 ,..., /i} such that the functions and af^ 
defined in (C.22) represent the limits e —)• 0 of the coordinates Og and ke defined via (C.15), 
respectively. 

The coordinates a = (oi,..., ah) may be completed into a system of Darboux coordinates 
{a,t) for AiniC) by introducing the coordinates t = {ti,... ,th) using a variant of the Abel 
map defined as 

^ rui 

tk = -Y. (C.23) 

1=1 

where Uh, k = 1,..., /i are the Abelian differentials of the first kind on the spectral curve 
S which are dual to the differentials Oj in the sense that = ^ik- The functions tr 
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represent coordinates on the Prym variety. The fact that the coordinates {a,t) represent 
Darboux coordinates for A4ii{C) follows from the fact that 

>V’"(a, u, z) 


d 

E 

1=1 ■ 


rui 


(C.24) 


is a generating function for the change of coordinates {u,v) ■(-)> (a, t). Indeed, note that 


1 d 

Wfc := — 

2,71 OCik 


(C.25) 


is an abelian differential on S satisfying f^.ojk = Sik as a consequence of (C.22). We may 
therefore conclude that W^{a,u, z) satisfies 


\ d d 

——W^{a,u,z) =tk, —W^{a,u,z) = -Vk, 
27r oak ouk 


(C.26) 

identifying W^(a, u, z) as the generating function for the change of coordinates (u, v) (a, t). 


D. Classical limits of conformal field theory 


We had in the main text introduced chiral partition functions u, r; b) and x, r; k) 

in Liouville theory and the WZWN model respectively. It will be helpful to parametrize the 
representation labels /3 and j appearing in the arguments of the functions Z^{l3,u,T;b) and 
Z'^^{j, X, r; k) as 

b^ = -, (D.I) 

£2 

k + 2 = (D.2) 


Q _ Q : ■ 

* ■ 2 

1 Oe 


Using this parametrization allows us to introduce chiral partition functions Z^(a, u, r; ei, € 2 ) 
and Z^'^{a, x, r; ei, 62 ) depending on two parameters ei and 62 - We may therefore define two 
different classical limits of Liouville theory and the S'L(2)-WZW model by sending ei or 62 
to zero, respectively. We are interested in the limit where both ei and 62 are sent to zero, but 
it helps to first study the limit ei —>• 0 with €2 finite before sending 62 —)• 0. After sending ei 
to zero we will find a relation to the moduli space of e 2 -connections. 

The two cases related to Virasoro and Kac-Moody algebra, respectively, can be treated 
in very similar ways. In each of these cases we will show that the leading asymptotic behavior 
of the chiral partition functions. 


logZ^^(a,x,r;ei,€ 2 )-—x, r; 62 ), 

£1 

log2:^(o,u,r;ei,e2) ~- y^{a,u,T;e2) 


(D.3) 
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is represented by functions x, r; € 2 ) and y^{a, u, r; € 2 ), which are generating functions 

for the changes of Darboux variables {x,p) <-)• (a,/c) and {u,v) -H- {a,k) for respec¬ 

tively. 

The dependence on the variables x (resp. u) will be controlled by the partial differen¬ 
tial equations satished by x, r; ei, 62 ) (resp. Z^{a,u,T]ei,€ 2 )), known as Knizhnik- 

Zamolodchikov-Bernard (KZB) and Belavin-Polyakov-Zamolodchikov (BPZ) equations. In 
order to control the dependence on the variables a in both cases the crucial tool will be the 
Verlinde loop operators defined by integrating the parallel transport defined by KZB- and 
BPZ-equations, respectively. The Verlinde loop operators can be represented as difference op¬ 
erators acting on the a-variables. The limit ei —)■ 0 of the relations between parallel transport 
and the corresponding difference operators will govern the a-dependence of y'^^{a, x, r; € 2 ) and 
y^{a, u, r; 62 ). The following discussion considerably refines the previous observations [67,68] 
by supplementing the “other side of the coin” represented by the Verlinde loop operators. 

To simplify the exposition we will spell out the relevant arguments only in the case when 
C has genus zero. The dependence on the complex structure of C may then be described using 
the positions z = (zi,..., Zn) of the marked points. We will therefore replace the parameters 
r by the variables z in the following. The generalization of this analysis to higher genus 
Riemann surfaces will not be too hard. 

D.l Preparations: Insertions of degenerate fields 

It will be useful to modify the conformal blocks by inserting a variable number of m extra 
degenerate helds at position y = (yi,... ,ym)- 

WZW model 

We will consider conformal blocks of the form 

y-X, z) := {^^'^{Xn\Zn) ... {xi\zi) ^^WmlVni) • • • {wi\yi))^^^. (D.4) 

We will impose the “null vector decoupling” equation on the degenerate field ^ 2 (^w\y): 

dl<^'2{w\y) = 0, (D.5) 

which means that transforms in the two-dimensional representation ~ C[rc]/(rc^) 

of SI 2 . It follows that Z^j^z(w,y',x, z) defines an element y'^^{y;x, z) of 

The corresponding chiral partition functions x, z) satisfy additional first order dif¬ 

ferential equations governing the y-dependence which will be formulated explicitly below. The 
family of chiral partition functions obtained in this way represents a convenient repackaging of 
the information contained in the chiral partition function Z^'^{x, z) without extra degenerate 
fields {m = 0). The chiral partition functions Z'^^{x,z) essentially represent the bound¬ 
ary conditions for the integration of the differential equations governing the y-dependence of 
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z). One may recover Z^^{x,z) from the family of x, z) by taking suitable 

limits. The presence of extra degenerate fields modifies the KZ-equations as 


ei OZr 


E 

r'=l 

r'jLr 


'Haa' 


rra rja 
Zy" Zy‘ ^ 


^'^^{y,X,z) + J2Vaa' 

S = 1 


'Taj.a 

Uy- 

Vs 


^'^\y;x,z), (D. 6 ) 


where denote the matrices representing SI 2 on the s-th tensor factor of and are 

the differential operators introduced in (B.3). In addition we get the following m differential 
equations: 

j.a^a' 

^^^{y;x,z) + y2 W ^ ^'^^{y;x,z), (D.7) 

s'^s 


ei oy. 


= E 

r=l 


Vaa' 


j.a >ja 
ig %Jy- 

Vs - Zr 


The space of solutions to the equations (D.7) is determined by the space of conformal blocks 
without extra degenerate fields m = 0. This follows from the fact that one may regard the 
partition function Z^^{x, z) as initial values for the solution of (D.7). One may, on the other 
hand, recover the partition functions Z^^{x,z) by considering even m and taking a limit 
where the insertion points ys collide pairwise. 


Liouville theory 

The situation is similar in the case of Liouville theory. In the presence of m degenerate fields 
of weight — ^ — 16 “^ and I degenerate fields of weight — ^ the chiral partition functions 

will satisfy I BPZ equations (D.9a) We shall consider the Liouville conformal blocks 


Z^{y,u,z) = ( n 


m I 

^OLr4>{Zr) n 




(D. 8 ) 


r=l 


s=l 


k=l 


C,ip 


The conformal blocks (D. 8 ) satisfy the null vector decoupling equations 




Ar 


+ 


1 


d 


i V (^/c ■2'r)^ '^k dz^ 

r=l ^ ' 


362 + 2 


d 


V4(rtfc - rifc/)2 Uk-Uk'duk' 


E 


k'=l 

fc+fc 


j^^\A{uk-ysY Uk-Vsdys 
\Zj^{y-yz) = 0 


(D.9a) 


1 ^2 
6 ^ dy 


+ E 


Ar 


+ 


1 


d 


E 

s'=l 

s'^s 


V {Vs - ZrY ys - Zr dz, 

352 + 2 1 d 


E 


36-2 2 


1 


d 


^{vs - Vs'Y Vs - Vs' dus' 


V4(ys - ys-Ukduk 

Z^{y,u,z) = 0. 


Equations (D.9) imply the fusion rules 

[V-h/ 2 \ • [D„] ~ [D„_+ 2 ] + [Va-b/ 2 ], 
[^-1/26] • [47.] ~ [V^-l/2b] + [47._/2b] . 


(D.9b) 


(D.IO) 

(D.ll) 
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D.2 Limit ei —)• 0 

We will next discuss the behavior of the solutions to the null vector decoupling equations in 
the limit ei —?■ 0 . 


WZW-model 

In order to study the limit ei —>• 0 it is useful to multiply (D.7) by ei and (D. 6 ) by ef. One 
may solve the system of equation (D. 6 ) and (D.7) with the following ansatz, 

n 

z){l + 0{ei)) , (D.12) 

S=1 


which will yield a solution to (D.7) provided ip{y]x,z) and y'^^{x,z) satisfy the following 
system of equations: 


d 

£2 A{y)]^l;{y;x,z) = 0 , 


where 




'Haa' 


fa A a' 

Lg 


Ar — 


r=l S'-- 


X'l'P^ ‘^L’pXj' x^Pj' 

Pf / t * X’pP’p 


(D.13a) 

(D.13b) 

(D.13c) 


We recognize model (B) for the flat connections. The limit of (D. 6 ) yields in addition 


d 


Hr :=e2-^y'^\x,z) = 


Vaa' 


r'=l 

r'jkr 


Aa Aa' 
Zn^ Zr-^ 


(D.13d) 


These equations characterize the Hamiltonians of the Schlesinger system. We have thereby 
reproduced results of [67,68]. 


Liouville theory 

In order to study the limit ei —)• 0 it is useful to multiply (D.9b) and (D.9a) by € 162 . One 
may solve the system of equation (D.9a) and (D.9b) with the following ansatz, 

n 

y^{y]u,z) = e~^^ JjT/>^(y;M,2;)(l + C>(ei)) , (D.14) 

S =1 

which will yield a solution (D.9b) provided u, z) and y^{u, z) satisfy the following system 
of equations: 


52 

+ tiys) )V’^(y; u,z)=0, 


(D.15a) 
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where 



Vk = -^^y^{u,z) , 
duk 


Sf — ; 


Vs 


Hr 



4(ys - Uk)"^ 


3e2 


Vs 



(D.lSb) 


(D.15c) 


The equations (D.9a) yield in addition 


OO 


vl + 4,2 = 0 > Kv) = ^ tk,i{y - UkY ^ , 


(D.16a) 



(D.16b) 


These equations define the Hamiltonians of the Gamier system. 

D.3 Verlinde loop operators 

The dependence of the chiral partition function on the variables a is controlled by the Verlinde 
loop operators. They are defined by modifying a conformal block by inserting the vacuum 
representation in the form of a pair of degenerate fields, calculating the monodromy of one 
of them along a closed curve 7 on C, and projecting back to the vacunm representation, 
see [9,69] for more details. A generating set is identified using pants decompositions. 

The calculation of the Verlinde loop operators is almost a straightforward extension 
of what has been done in the literature. The necessary results have been obtained in [9, 
69] for Liouville theory without extra insertions of degenerate fields H_b/ 2 (y)- It would be 
straightforward to generalize these observations to the cases of our interest. For the case 
of Kac-Moody conformal blocks one could assemble the results from the known fusion and 
braiding matrices of an extra degenerate held ^ 2 (^iu,y). As a shortcut let us note, however, 
that the results relevant for the problem of our interest, the limit ei —)• 0 , can be obtained in 
a simpler way. 

One may start on the Liouville side. The key observation to be made is the fact that the 
presence of extra degenerate helds V_i/ 2 b{y) modihes the monodromies of V_f,/ 2 {y) only by 
overall signs, as the monodromy of V_fe/ 2 (y) around V_ii 2 b{uk) is equal to minus the identity. 
It is nseful to observe (see Appendix E.4) that the separation of variables transformation maps 
the degenerate held $2 (tf;^ y) to the degenerate held V_f,/ 2 {y)- It follows that the monodromies 
of ^ 2 (^w,y) must coincide with the monodromies of V_i,/ 2 {y) up to signs. Using the results 
of [9,69] we conclude that 



(D.17) 
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where £ {=tl} ^ind £ {iljj while the explicit expressions for the difference operators 
Le,s; nre 


-e,s = 2cosh(27rae/e2). 


Le,t = 


2 COs(7rei/€2)(Le^2Le,3 “1“ Le^lLe^4) -|- (Lg^l+ Lg 2Lg,4) 
(ag + iei))2sinh(f^(ag-iei)) 

_ p-Kike/e2 \/ci2(Lr,g)c34(Lr,g) ^^fk./e-;> 

2sinh(27rae/e2) 


2sinh(^ 
'' £2 

1 


(D.18a) 

(D.18b) 


+ V _ 

v/2 sinh(27rae/e2) 


^2 sinh(27rae/e2) ' 


using the notation Cij{\-e,s) 


\2 + l2 . + . + L 

e,s ^ e,* ^ ej ^ e 


-4, and 


€162 Q 

27ri dae 


(D.19) 


As the KZB-equations (D.7) turn into the horizontality condition (D.13a), the Verlinde 
loop operators will turn into trace functions when ei —)■ 0. The limit of the left hand side of 
(D.17) is therefore found by replacing 7r''(7g^s) and 7r^(je,t) with the expressions in (C.15), 
calculated from the connection A(i/) appearing in (D.13a). Note that the connection A(y) is 
thereby defined as a function of the parameters x and a. The limit ei —>• 0 of the right hand 
side of (D.17) is straightforward to analyze by using (D.3) and (D.18). It can be expressed 
in terms of the derivative of with respect to the variable a. In this way one finds that 
the the limit ei —)■ 0 of equations (D.17) implies the relations 

ke{a,u) = e 2 -^^y^'^{a,u,z). (D.20) 

2ilT 

Equation (D.20) identifies y^'^{a, u, z) as the generating function for the change of variables 
{x,p) {a,k). The analysis in the Liouville case is very similar. 


D.4 Limit e 2 — )■ 0 

It remains to discuss the behavior in the limit €2 —>• 0 of y'^^(a, x, z; € 2 ) and y^{a,u, z; € 2 )- 
We claim that in the two cases we find a behavior of the form 

y^'^{a,x,z) ~ —y^^{a,z)+W^^{a,x,z) + ..., 

£2 

y^{a,u,z) ~ —y^(a,z)+W^{a,u,z) + ... , 

£2 

where F^'^{a,z) = F^{a,z), while W^^{a, x, z) and W^(a, x, z) are the generating functions 
for the changes of variables {x,p) •(-)• (a, t) and {u,v) -H- (a, t), respectively. 

We begin by considering (D.22). The equation (D.15) can be solved to leading order by 
a WKB-ansatz 

k(yw,^) - {v{y)f + t{y) = 0. (D.23) 


(D.21) 

(D.22) 
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The asymptotics of the generating function y^{a, u\ z) which coincides with the classical 
Liouville conformal blocks will be of the form 

1 


y^{a, u, z) ~ z) + >V^(a, u,z) + ... , . 

€2 

Indeed, an expansion of the form will satisfy (D.16) and (D.20) if J-i^{a,z) satisfies 

d 




D 

e ’ 


(D.24) 


(D.25) 


identifying y(a,z) as the prepotential, and if furthermore 


d 

duk 


>V^(a, u, 


z) = -Vk- 


(D.26) 


This means that 


>V^(a, u; z) 



V . 


Following the discussion in Appendix C.4 we may identify yV^{a,u] z) as the generating 
function of the standard change of Darboux variables (u, v) -H- (a, t) which is defined by the 
Abel map. 

The corresponding statement for yV'^^{a,x; z) now follows easily from (D.13c), and the 
fact that y^^{a, X, z) and y^{a,u,z) differ only by the generating function [x] u, z) for 
the change of Darboux variables {x,p) -H- {u,v) which does not depend on a. 


E. Explicit relation between Kac-Moody and Virasoro conformal blocks 

We will explain in this appendix how to obtain an explicit integral transformation between 
the conformal blocks in Liouville theory and in the WZW model using the observations made 
in Section 4.7. This is the separation of variables (SOV) relation (1.1) which we discussed in 
the Introduction. 

E.l SOV transformation for conformal blocks 

In order to partially fix the global sb-constraints we shall send Zn ^ oo and Xn oo, 
defining the reduced conformal blocks Z^^{x,z) which depend on x = (xi,..., x„_i) and 
z = (zi,..., Zn-i). Let Z'^^{p,z) be the Fourier-transformation of the reduced conformal 
block Z'^^{x, z) of the WZW model w.r.t. the variables x. It depends on = {pi, ..., pn-i) 
subject to = 0. There then exists a solution Z^{y, z) to the BPZ-equations 

Vll^-Z^ = 0, VA: = 1,...,/, (E.l) 

with differential operators given as 

pBPz ^ 52 ^ ^ A _ f + 2 _ 1 (9 \ 

9ul ^^KiUk-Zr)"^ Uk-ZrdZrJ ^^yAiuk-Uk')"^ Uk - Uk' dUk') ' 

k'^k 
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such that the following relation holds 

= uq 5{YIIZI fJ'i) Qn{y,z)Z^{y,z) 
The function Qniy,z) that appears in this relation is defined as 


(E.2) 


Qniy,z)= {Zr - Zs)^ 

r<s<n—l k<l<n—3 


n—1n—3 


1 


{Uk - Ul)^ Y\ YlZr - Uk) 

r=l k=l 

The relation (E.2) will hold provided that the respective variables are related as follows: 
(1) The variables yi ,..., yn-i are related to tti,..., Un-S; via 


(E.3) 


n—1 

E 

r=l 


fJjr 


nizlit-uk) 


t — Zr 


= Uq - 


n:ziit-zr) 


(E.4) 


In particular, since YZ^=i = 0, we have uq = Z2r=i HrZr- 

( 2 ) b^ = -{k + 2)-\ 

(3) The Liouville momenta are given by 

ar = a{jr) ■■= b{jr + 1 ) + 


(E.5) 


We may use formula (E.2) to construct bases of solution to the KZ-equations from Liouville 
conformal blocks. 

E.2 Reformulation as integral transformation 

We want to write the expression for Z'^^(x,z) 


n—1 


Z^^(x,z) = f ^ 6(J2r^i Mr) &n(u,x) Zi^(u,z) Yl 


Jr 


r=l 


Ml Mn—l 

as explicitly as possible. To this aim let us note first that 

/ \ \ l \ \ / \ Ofc=l 

Mr{u) = UoXr[u) , Xr{u) := - 


m^Zr-Zs) ’ 


and furthermore 

d/il dfin—l r-z-^n-l \ r\ t \ \ ^ 

- . . . - '^(Lr=l Mr) &n{u\x) = - du{u) 

Ml Mn-1 Uq 

n—1 n—1n—3 

(Zr - 


(E.6) 


(E.7) 


(E.8) 


n—3 


k<l 


dv{u) := dui. .. dUn-3 (Zr - Zs)^'^^ n JJ (^r - 'Wfc) ^ ^11 “ Ul) 

r^s r=l k=l 

n—1 ^ n—3 

dui...dUn-3 ]^[Ar(u)]“^"5b^ - tt;) 


1 +; 


\ 1 +; 


r=l 


k<l 
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We may therefore calculate 


n—1 r. 1 n—1 

^^0 —J ^ I I ^IUqXtXt 


= [ dv{u) ITA^-^’" [ Uq-^ Z^{u,z)\\> 

J .=1 J ^0 .=1 


where J := —jn + YXi=i Jr- The integral over uq is of the form 


7 n—L 

,,-J JuoXrXr _ 


uo 


N.,{Y.>' 


r=l \ r=l 

where Nj depends neither on x nor on z. It follows that 

J 


.j'tJU'f' 1 ^ 


n—1 


\x,z) =Nj du{u) (^X]r=l ArXr^ Zl(u, 2 ;)]^A, 


~3r 
r 5 


r=l 


= Nj / dui...dUn-3 w) z) , 


where the kernel IC^°'^{x,u) is defined as 
}C^°^{x,u) := 


(E.9) 


(E.IO) 


(E.ll) 


(E.12) 


n-l T-rn-S 


lXk=li^r - Uk) 

/ V ^ y-rn—l/ \ 

[[s^r\Zr-Zs) 


n—3 


n—1 


ii{uk-uiX+^ n 

k<l r=l L 


ngk^r-^.) 

nj;f(2, - u,) 


GLrjh 


Note that the x-dependence it entirely in the first factor on the right hand side of (E.12). 

The choice of contours in (E.ll) is a delicate issue that we will not address here. Using 
the standard contour M in the definition of the Fourier-transformations in (E.6) will of course 
determine a particular choice of contours in (E.ll). Any choice of contours that ensures 
absence of boundary terms in the relation between the differential equations satisfied by 
Z^'^ix^z') and Z^iu^z') could also be taken to define a relation of the form (E.ll) between 
bases of conformal blocks in the WZW-model and in Liouville theory. Changing the contours 
in (E.ll) amounts to a change of basis in the space of solutions to the KZ-equations obtained 
from a fixed basis in the space of Liouville conformal blocks. It would be interesting to 
identify the basis defined by (E.ll) for a given choice of contours precisely, and to investigate 
the dependence on the choice of contours. 


E.3 Semiclassical limit 

Now we consider semiclassical limit ei, 62 —>■ 0, setting 

Oir = (eie2)"2 . 


We have then 

loglC^°^ix,u) = er^>V®°^(x,u) + 0(e?), 


(E.13) 

(E.14) 
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with 


li) = Klog 


n-l T-rn-a 






r 

n-l 


(E.15) 


+ '^lr 


n—1 


n—3 


y~l ^Og{Zr - Zs) log{Zr - Uk) 


k=l 


r=l L s^r 

We have denoted k := —In + X]r=i If we send only ei —>• 0, we get a modified result: 

logfC^°^ix,u) = er^>V®°''(x,M;e 2 ) + 0(e?) (E.16) 


with 


W®°'^(x,n;e2) = ^log 


n-l T-rn-S 


nUi^r -Uk) 


/ V ■■' T—rn—i / 

1^1 Us^ri^r-Z. 


+ jY^^ogiuk - Ul) 


n—3 


(E.17) 


k<l 


n—1 


+ 


r=l L s^r 


n—1 


n—3 


log(2;r - ^s) - log(Zr - Mfc) 


fc=l 


E.4 SOV transformation in the presence of degenerate fields 

We also use the version of this correspondence in the presence of the fields <h 2 (u)|y), as 
appear in (D.4). This is kind of interesting. Note that the Eourier-transformation of the null 
vector equations d^ 4 > 2 (rt;|y) = 0 gives //^<h 2 (/i|y) = 0. This indicates that conformal blocks 
containing <f> 2 (^|y) must be understood as distributions with support at /r = 0. If we send 
/Tr —^ 0 in the change of variables (E.4), we will loose the pole at t = Zr on the left hand side. 
This means that one Uk must approach Zr in order to cancel the pole at t = Zr on the right 
hand side of (E.4). It follows that the degenerate field e~^ V(“fc) fuses with the field _ 

Applying these observations to the case where jr = 1/2, which corresponds to = \Q + h we 
get as leading term in the OPE of e~^ <t>iuk)^2ar<j>{zr) ^ ffeld with conformal dimension A_fe/ 2 ) 
which is degenerate. This indicates that the WZW conformal blocks (D.4) can be represented 
in terms of the Liouville conformal blocks (D.8) with I = n — 3, where the insertion of a field 
4 > 2 (t(;^|y^) corresponds to the insertion of Even if the argument above may look 

delicate, the conclusion seems hard to avoid: We need to map the field $ 2 (t(;|y) to another 
field with two-dimensional monodromy. The only candidate with the right behavior for 6 —)• 0 

is 
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